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Abstract 

Adaptive  optics  applies  advanced  sensing  and  control  to  improve  the  ability  of 
optical  systems  to  collect  images  through  a  turbulent  atmosphere.  The  results  of  this 
research  effort  demonstrate  that  the  combination  of  two  recent  approaches  improves 
the  performance  of  adaptive  optics  in  directed  energy  and  laser  communication  sce¬ 
narios.  The  first  approach  is  adaptive  control,  which  offers  improved  performance 
over  fixed-gain  controllers  in  the  presence  of  rapidly  changing  turbulence.  The  second 
approach  incorporated  into  the  study  is  a  dual-mirror  system.  The  two  mirrors  are 
a  high-bandwidth,  low-actuator-stroke  (tweeter)  mirror  and  a  low-bandwidth,  large- 
actuator-stroke  (woofer)  mirror.  The  woofer-tweeter  combination  allows  for  better 
compensation  of  the  large- variance,  high-spatial-frequency  phase  distortion  generated 
by  strong  turbulence.  Two  different  adaptive  controllers  are  presented,  one  using  a 
relatively  simple  model  reference  adaptive  system  controller  and  one  using  a  lattice 
filter  controller.  The  lattice  filter  is  implemented  in  two  ways.  In  one  implemen¬ 
tation  the  filter  operates  on  the  individual  actuators,  while  in  the  other  it  operates 
on  frequency  weighted  modes.  The  modal  implementation  reduces  the  computational 
burden  of  the  filter.  The  performance  of  the  different  adaptive  controllers  is  compared 
to  both  each  other  and  to  a  traditional  fixed-gain  controller.  Simulations  show  that 
adaptive  control  of  woofer-tweeter  AO  can  increase  the  mean  Strehl  ratio  by  up  to 
20%.  In  general,  the  lattice  filter  controllers  outperform  the  model  reference  adaptive 
system  controller.  However,  in  cases  where  the  lattice  filter  cannot  use  a  sufficient 
number  of  modes,  the  model  reference  adaptive  system  can  outperform  the  lattice 
filter. 
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Adaptive  Control  Of 


Woofer-Tweeter  Adaptive  Optics 

I.  Introduction 

Adaptive  optics  (AO)  refers  to  a  closed-loop  optical  system  whose  purpose  is 
to  correct  for  atmospheric  aberrations.  In  astronomical  applications,  adaptive 
optics  is  used  to  reduce  the  distortion  of  the  incoming  light  and  improve  the  image. 
In  laser  communication  and  directed-energy  applications,  an  AO  system  pre-distorts 
the  outgoing  light.  Then,  the  atmosphere  corrects  the  pre-distortion  which  results  in 
a  more  focused  spot  at  the  target  and,  consequently,  more  energy  on  the  target. 

There  are  two  recent  approaches  to  improve  the  performance  of  adaptive  op¬ 
tics  in  directed  energy  and  laser  communication  scenarios.  The  first  approach  is 
adaptive  control,  which  offers  improved  performance  over  fixed-gain  controllers  in  the 
presence  of  rapidly  changing  turbulence.  The  second  approach  is  a  dual-mirror  sys¬ 
tem.  The  two  mirrors  are  a  high-bandwidth,  low-actuator-stroke  (tweeter)  mirror  and 
a  low-bandwidth,  large-actuator-stroke  (woofer)  mirror.  The  woofer-tweeter  combi¬ 
nation  allows  for  better  compensation  of  the  large- variance,  high-spatial-frequency 
phase  distortion  generated  by  strong  turbulence.  This  research  combines  these  two 
approaches. 

1.1  Problem  Statement 

The  objective  of  this  research  is  to  develop  an  adaptive  controller  for  a  woofer- 
tweeter  AO  system  which  improves  the  performance  of  the  system  compared  to  tra¬ 
ditional  fixed-gain  controllers.  The  controller  design  is  verified  using  wave-optics 
simulations  in  different  atmospheric  turbulence  scenarios. 
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The  particular  goals  of  this  research  are: 


•  to  determine  a  control  architecture  suitable  for  implementing  an  adaptive  con¬ 
troller  with  a  woofer-tweeter  AO  system, 

•  to  develop  an  adaptive  controller  using  the  determined  control  architecture, 
which  improves  the  performance  of  a  woofer-tweeter  AO  system, 

•  to  reduce  the  computational  complexity  of  the  controller  through  the  use  of 
modal  decomposition, 

•  to  use  wave-optics  simulations  to  demonstrate  the  improvement  gained  by  using 
an  adaptive  controller, 

•  and  to  compare  the  improvement  with  various  implementations  of  the  adaptive 
controller. 

This  research  effort  realized  all  of  the  stated  goals.  The  analysis  shows  that 
an  adaptive  controller  improves  the  performance  of  woofer-tweeter  AO  systems  and 
adaptive  control  of  woofer-tweeter  AO  is  a  viable  approach  to  compensating  strong, 
fast-moving  turbulence. 

1.2  Thesis  Overview 

Chapter  II  discusses  the  topic  of  atmospheric  turbulence  and  its  mathemati¬ 
cal  description.  It  then  provides  a  brief  overview  of  AO  systems  and  their  major 
components,  including  commonly  used  controllers.  After  covering  typical  AO  sys¬ 
tems,  woofer-tweeter  AO  systems  are  discussed.  First,  the  need  for  woofer-tweeter 
AO  is  established  followed  by  a  discussion  on  woofer-tweeter  controllers.  The  last 
section  in  Chapter  II  defines  adaptive  control,  covers  two  types  of  adaptive  control 
techniques,  and  discusses  current  applications  of  adaptive  control  and  AO,  including 
modal  decomposition. 
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Chapter  III  gives  details  on  the  parameters  used  to  model  the  atmospheric  tur¬ 
bulence  and  the  AO  system  hardware.  These  parameters  are  needed  to  derive  the 
specific  control  laws  used  in  the  adaptive  controllers,  as  well  as  define  the  simula¬ 
tions.  With  the  required  system  parameters  established,  three  different  controllers 
are  derived  using  the  techniques  discussed  in  Chapter  II. 

Chapter  IV  implements  the  controllers  designed  in  Chapter  III.  First,  a  baseline 
for  comparing  the  different  controllers  is  established  using  a  standard  woofer-tweeter 
system.  Each  controller  is  then  simulated.  The  simulations  are  used  to  determine  the 
best  implementation  for  each  controller.  Once  the  specific  implementation  has  been 
established,  the  controllers  are  compared  to  the  baseline  system  and  to  each  other. 

Chapter  V  summarizes  the  simulation  results  and  provides  a  detailed  comparison 
between  the  different  controllers.  Based  on  the  evidence  presented,  a  final  conclusion 
is  given  that  adaptive  control  indeed  improves  woofer-tweeter  AO.  Various  areas  for 
future  work  are  then  recommended.  These  recommendations  propose  ways  to  use  or 
build  upon  the  adaptive  control  algorithms  presented. 
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II.  Background  and  Current  Research 


This  chapter  gives  a  brief  overview  of  conventional  adaptive  optics  systems,  in¬ 
cluding  the  major  components  and  common  control  designs.  Two  areas  of  recent 
interest  to  the  AO  community,  woofer-tweeter  AO  and  adaptive  control,  are  also  dis¬ 
cussed.  The  topics  presented  in  this  chapter  form  the  foundation  of  the  adaptive 
woofer-tweeter  controller  design  presented  in  Chapter  III. 

2.1  Conventional  AO 

A  typical  AO  system  uses  the  following  components:  a  wavefront  sensor  (WFS) 
to  measure  the  light  from  a  source  beacon,  a  deformable  mirror  (DM)  to  adjust 
the  waves,  and  a  control  computer  to  translate  the  WFS  measurements  into  DM 
commands.  To  better  understand  AO,  the  following  four  sections  discuss  atmospheric 
turbulence,  deformable  mirrors,  wavefront  sensors,  and  control.  Figure  2.1  shows  how 
these  components  are  used  in  a  typical  AO  system. 

2.1.1  Atmospheric  Turbulence.  In  order  to  understand  AO  systems,  it  is 
necessary  to  understand  how  the  atmosphere  affects  light  and  how  to  describe  and 
quantify  these  effects.  Because  air  has  an  index  of  refraction  that  is  not  exactly  one, 
light  is  bent  as  it  travels  through  the  air.  The  index  of  refraction  depends  on  the 
pressure,  density,  and  temperature  of  the  air  as  well  as  the  frequency  of  the  light.  In 
a  given  mass  of  air,  the  pressure,  density,  and  temperature  vary  throughout  the  mass. 
These  variations  result  in  pockets  of  differing  indices  of  refraction.  The  covariance 
of  the  refractive  index,  Bn ,  can  be  computed  by  comparing  the  refractive  index  at 
one  point,  7ii(rq),  with  another  point,  n\(r  +  rq),  and  averaging  over  the  air  mass  of 
interest  according  to 

Bn{r)  =  (: m(r  +  rq)  ni(rq))  (2.1) 

where  rq  is  the  position  in  the  air  mass,  r  is  the  distance  between  the  points,  and  (•) 
is  the  expectation  operator.  The  power  spectral  density,  $(iF),  of  the  atmosphere  is 
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Incident  Wavefront 


Figure  2.1:  A  typical  AO  system  including:  a  wavefront  sensor  (WFS),  deformable 
mirror  (DM),  and  a  control  computer  [21]. 

the  three-dimensional  Fourier  transform  of  the  covariance  give  by 

$(K)  =  Ji  Bn(r )  exp  (—A"  ■  r)d3r  (2.2) 

where  l0  is  the  inner  scale,  or  smallest  atmospheric  eddy  size,  A0  is  the  outer  scale,  or 
largest  atmospheric  eddy  size,  and  K  is  the  three-dimensional  spatial-frequency  [20]. 
In  general,  this  integral  does  not  have  a  closed-form  solution.  However,  with  certain 
assumptions  a  closed-form  solution  can  be  found.  Two  of  the  most  common  solutions 
are  the  Kolmogorov  spectrum  and  the  Von  Karrnan  spectrum.  In  the  Kolmogorov 
spectrum,  the  inner  scale  is  zero  and  the  outer  scale  is  finite,  which  yields  the  following 
expression  for  the  power  spectral  density: 

$(AT)  =  0.033C^A'-11/3  (2.3) 
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In  the  Von  Karman  spectrum,  the  inner  scale  is  zero  and  the  outer  scale  is  finite  and 
yields  the  following: 


$(K)  =  3.9  x  10"5C^Lj1/3 


LlK2 
1  +  — — 
(2tt)2 


-11/6 


(2.4) 


In  both  equations  the  atmospheric  structure  constant,  C2,  represents  the  severity 
of  the  turbulence  and  typically  ranges  between  10-15  m-2/3  and  10~18  m-2/3  as  a 
function  of  altitude  [20].  For  simulation  purposes,  several  canonical  models  of  C2 
exist  within  the  AO  community  [20]. 

The  various  pockets  of  air  with  different  refractive  indices  can  be  thought  of 
as  a  collection  of  lenses  floating  in  the  air.  Coherent  light  passing  through  these 
eddies  becomes  spatially  decorrelated  at  an  aperture  because  light  arriving  at  different 
points  of  the  aperture  travels  through  different  paths.  A  common  metric  describing 
the  strength  of  the  turbulence  is  Fried’s  parameter,  r0.  This  parameter  refers  to 
the  largest  aperture  over  which  light  would  be  coherent  for  a  given  turbulence.  In 
other  words,  a  benign  atmosphere  with  no  turbulence  would  have  r'o  =  oo.  As  the 
atmosphere  becomes  more  turbulent,  the  value  of  r0  decreases.  Thus,  smaller  values 
of  r0  indicate  a  more  severe  turbulence.  If  Kolmogorov  turbulence  is  assumed,  the 
coherence  length  for  a  plane  wave  is 


r0  =  1.68  (C^Lk2)  3/5  (2.5) 

where  L  is  the  length  of  the  propagation  path  and  k  is  the  wave  number.  The  wave 
number  is  related  to  the  wavelength,  A,  by 

k  =  f  (2.6) 
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In  astronomical  applications  and  other  scenarios  with  vertical  propagation,  the  struc¬ 
ture  constant,  C%,  varies  with  altitude,  z,  and  zenith  angle,  /3.  In  this  case, 


=  0.423/c2  sec (3  /  C2(z)dz 


for  a  plane-wave  source.  Typical  values  of  r0  range  from  5-10  cm. 

Another  metric  to  describe  atmospheric  turbulence  is  the  isoplanatic  angle,  90. 
The  isoplanatic  angle  is  the  angular  radius  around  the  beacon  that  exhibits  the  same 
turbulence  as  the  beacon.  In  other  words,  an  object  within  9q  of  the  beacon  will  be 
properly  imaged  by  the  AO  system.  At  a  given  zenith  angle,  90  is 


2.91k2  sec8/3  (3  /  C2{z)z5/3dz 


Typical  values  of  90  are  between  7-10  //rad. 

All  of  the  metrics  presented  above  describe  spatial  characteristics  of  turbulence. 
One  measure  of  the  temporal  characteristics  of  turbulence  is  the  Greenwood  frequency, 
fa-  The  inverse  of  the  Greenwood  frequency  is  Greenwood  time  constant,  To-  This 
time  constant  describes  the  interval  over  which  turbulence  is  nearly  constant  and  is 


given  by 


r0=  2.91k2  sec  (3  /  C2(z)V5/3(z)dz 


where  V (z)  is  the  wind  velocity  along  the  path.  If  the  wind  velocity  is  constant  then 


To  =  2.34^ 


(2.10) 


or  equivalently 


V 

fa  =  0.427- 
r  o 


(2.11) 


Typical  time  constants  are  on  the  order  of  milliseconds  [2], 
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The  performance  of  an  AO  system  can  be  measured  by  how  well  it  compensates 
for  the  atmospheric  turbulence.  A  perfectly  compensated  wavefront  would  be  a  plane 
wave  with  uniform  phase  across  the  aperture.  If  this  perfect  plane  wave  were  imaged 
it  would  result  in  a  spot  whose  intensity  is  determined  by  the  optical  properties  of 
the  imaging  system.  This  is  known  as  the  diffraction-limited  case.  An  aberrated 
wavefront  results  in  a  less  focused  spot  with  reduced  intensity.  The  ratio  of  the  peak 
aberrated  spot  intensity  to  the  peak  diffraction-limited  intensity  is  known  as  the  Strehl 
ratio.  The  Strehl  ratio  is  always  between  zero  and  one  and  is  used  to  measure  an  AO 
system’s  performance.  An  AO  system  which  perfectly  compensates  the  atmospheric 
distortion  would  have  Strehl  ratio  of  one.  As  the  wavefront  becomes  more  distorted, 
the  Strehl  ratio  approaches  zero. 

Using  the  atmospheric  turbulence  models  and  the  metrics  described  in  this  sec¬ 
tion,  it  is  possible  to  build  computer  simulations  with  the  appropriate  characteristics 
and  measure  the  performance  of  a  simulated  AO  system.  The  simulation  scenar¬ 
ios  presented  in  Section  3.1.2  are  described  in  terms  of  the  Fried  parameter  and  the 
Greenwood  frequency.  The  characteristics  of  the  atmosphere  are  also  essential  to  un¬ 
derstanding  the  need  for  two  DMs  as  discussed  in  Section  2.2.  In  order  to  understand 
how  an  AO  system  compensates  for  turbulence,  the  following  sections  describe  the 
components  of  an  AO  system. 

2.1.2  Deformable  Mirrors.  The  deformable  mirror  in  an  AO  system  allows 
the  system  to  alter  the  phase  of  the  light  in  response  to  the  WFS  measurements.  The 
goal  of  the  AO  system  is  to  adjust  the  figure  of  the  DM  so  it  has  the  same  shape 
as  the  aberrated  wave  with  the  opposite  sign.  If  the  aberrated  wave  is  of  the  form 
Aexp(i0),  then  the  DM  ideally  imparts  the  complex  conjugate,  Aexp(— if),  to  cancel 
the  aberration. 

A  deformable  mirror  consists  of  a  reflective  faceplate  with  actuators  below  the 
mirror  surface  to  move  it  up  or  down.  The  actuators  change  the  path  length  trans- 
versed  by  the  light,  which  changes  the  phase.  Because  the  light  reflects  off  the  mirror, 


the  change  in  the  total  distance  traveled  by  the  light  is  twice  the  displacement  of  the 
actuator,  once  before  reflection  and  once  after  reflection.  This  change  in  the  distance 
the  light  travels  is  known  as  the  optical  path  difference  (OPD).  Most  DMs  fall  into 
one  of  two  categories:  segmented  and  continuous  faceplate.  For  a  DM  with  a  seg¬ 
mented  faceplate,  the  surface  is  not  continuous  and  consists  of  many  small  mirrors 
placed  close  together.  Each  small  mirror  sits  atop  an  actuator  which  moves  it  up  and 
down.  Due  to  the  segmented  design,  each  actuator  influences  only  the  position  of  the 
corresponding  segment. 

In  a  continuous  faceplate  DM,  the  mirror  surface  is  a  continuous  sheet  with 
an  array  of  actuators  below  it.  As  a  result,  when  an  actuator  moves,  its  influence 
is  not  limited  to  the  mirror  region  directly  above  it.  The  spatial  deformation  of 
the  mirror  surface  due  to  an  actuator  movement  is  known  as  the  influence  function. 
The  amount  of  deformation  that  extends  to  an  adjacent  actuator  position  is  called 
coupling.  Figure  2.2  illustrates  the  differences  between  the  two  mirror  types. 


Segmented  Deformable  Mirror 


Figure  2.2:  There  are  two  types  of  DMs:  (a)  segmented  and  (b)  continuous  face¬ 
plate  [21]. 
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Due  to  their  lack  of  coupling,  segmented  mirrors  are  able  to  correct  for  aber¬ 
rations  with  high  spatial-frequency  content,  provided  the  individual  actuators  are 
sufficiently  close.  They  do,  however,  suffer  from  signal  loss  as  light  passes  through 
the  gaps  in  the  segments  and  from  diffraction  effects  at  their  edges.  A  continuous 
faceplate  DM  provides  a  smooth  fit  to  the  wavefront  as  a  result  of  coupling,  which  is 
advantageous  when  correcting  lower-order  aberrations. 

2.1.3  Wavefront  Sensors.  In  order  to  control  the  DM  and  correct  aberra¬ 
tions  in  phase  of  the  wavefront,  it  is  necessary  to  measure  the  light.  A  wavefront 
sensor  (WFS)  provides  this  measurement.  Ideally,  a  sensor  which  directly  measures 
the  phase  of  the  wavefront  would  be  used.  Unfortunately,  direct  measurements  of  the 
phase  are  not  possible.  Several  methods  exist  which  allow  the  phase  to  be  inferred 
from  other  measurements.  The  most  common  WFS,  a  Shack- Hartmann  WFS,  uses  a 
non-interferometric  approach  and  measures  the  spatial  derivative,  or  tilt,  of  the  phase. 
There  are  numerous  methods  which  use  various  implementations  of  interferometry  to 
create  diffraction  patterns.  The  intensity  of  the  diffraction  pattern  is  proportional  to 
the  phase  of  the  wavefront.  One  of  the  more  recent  implementations  is  known  as  the 
self-referencing  interferometer.  This  subsection  will  describe  the  important  elements 
of  these  two  types  of  wavefront  sensors. 

2. 1.3.1  Shack-Hartmann  WFS.  A  Shack- Hartmann  WFS  (SH  WFS) 
is  based  on  a  classical  optical  test  known  as  the  Hartmann  test.  In  a  Hartmann  test, 
the  wavefront  of  a  beam  is  measured  using  an  array  of  pinholes  to  determine  the 
slope  of  the  wavefront  over  a  number  of  points.  If  the  wave  were  perfectly  planar  each 
pinhole  would  produce  a  spot  directly  behind  it.  However,  if  there  are  aberrations 
in  the  wave,  they  cause  the  spot  position  to  be  shifted  by  an  amount  proportional 
to  the  wavefront  slope  at  the  pinhole  location.  Instead  of  an  array  of  pinholes,  a 
SH  WFS  uses  an  array  of  small  lenses,  called  lenslets,  to  focus  the  wavefront  into 
spots.  By  using  a  lenslet  array  which  captures  more  light,  the  SH  WFS  does  not 
suffer  the  same  signal  loss  as  a  pinhole  array  which  only  utilizes  the  small  fraction  of 
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Incident  Optical  Field 


Lenslet  Array 


Figure  2.3:  Example  of  a  SH  WFS  using  a  one-dimensional  lenslet  array.  The 

lenslet  array  produces  a  spot  on  the  detector  array.  The  spatial  shift  in  the  spot  is 
proportional  to  the  slope  of  the  wavefront  [21]. 

light  which  passes  through  the  pinholes.  Figure  2.3  shows  a  one-dimensional  example 
of  a  SH  WFS  where  each  lenslet  in  a  lenslet  array  focuses  a  portion  of  the  wavefront 
to  produce  a  spot  on  the  detector  plane.  Each  spot  is  displaced  from  the  center  of  the 
corresponding  section  of  the  detector  plane  due  to  the  local  slope  of  the  wavefront. 

One  method  for  determining  the  shift  in  the  spot  location  is  by  using  a  quadccll 
detector.  A  quadcell  consists  of  four  closely  spaced  detectors,  as  shown  in  Figure  2.4. 
The  spot’s  horizontal  displacement  is  calculated  by  averaging  spot  intensity  of  the 
two  left  detectors  and  subtracting  the  average  intensity  of  the  right  two  detectors. 
Similarly,  the  vertical  displacement  is  the  difference  between  the  averages  of  the  top 
and  bottom  pairs  of  detectors. 

Because  the  SH  WFS  does  not  provide  a  direct  measurement  of  phase,  it  is 
necessary  to  reconstruct  the  phase  from  the  WFS  measurements  in  order  to  drive 
the  DM.  To  perform  the  reconstruction,  it  is  essential  to  know  how  the  WFS  and 
DM  are  aligned.  The  alignment  geometry  refers  to  the  manner  in  which  the  WFS 
measurements  relate  to  the  DM  actuators.  Figure  2.5  shows  four  of  the  most  common 
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Quad  Cell 


Figure  2.4:  A  quadcell  detector  is  used  in  a  SH  WFS  to  measure  the  vertical  and 
horizontal  shift  of  the  focused  spot.  The  horizontal  displacement  is  the  difference  be¬ 
tween  the  left  and  right  average  intensities.  The  vertical  displacement  is  the  difference 
between  the  top  and  bottom  average  intensities  [21]. 

geometries.  The  Fried  geometry  measures  orthogonal  slopes  centered  between  four 
actuators.  The  Wavefront  Control  Experiment  (WCE)  geometry  is  identical  to  the 
Fried  geometry,  but  is  rotated  45°.  The  Hudgin  geometry  measures  orthogonal  slopes 
between  adjacent  actuators.  The  Southwell  geometry  measures  orthogonal  slopes  at 
each  actuator  position. 

Once  the  geometry  has  been  established,  the  relationship  between  measurement 
vector,  y,  and  the  actuator  commands,  c,  can  be  defined.  Assuming  a  linear  relation¬ 
ship  between  the  measured  wavefront  slopes  and  actuator  commands,  the  relationship 
is  simply 

Y(Mxl)  =  F(mxN)c(Nx1)  (2-12) 

where  the  matrix  F,  the  poke  matrix,  represents  the  system’s  geometry,  actuator  cou¬ 
pling,  and  influence  functions  and  describes  how  the  actuator  commands  correspond 
to  the  WFS  measurements.  The  subscripts  denote  the  dimensions  of  the  variables. 
In  general,  M  >  N  and  the  system  is  over  determined.  When  the  system  is  overde¬ 
termined  a  common  method  of  solving  Equation  (2.12)  is  to  use  the  least-squares,  or 
pseudo-inverse  solution, 

c=  (rTr)“1rTrj/  (2.13) 
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Figure  2.5:  Common  actuator- WFS  geometries:  (a)  Fried,  (b)  WCE,  (c)  Hudgin, 
(d)  Southwell.  The  numbered  circles  show  the  positions  of  four  DM  actuators,  the 
light  blue  circles  represent  the  subapertures  measured  by  the  sensor,  and  the  arrows 
show  how  the  slopes  are  defined  in  those  subapertures  [20]. 


13 


Because  the  relationship  between  phase  and  actuator  position  is  not  unique,  the 
matrix  (r2T)_1  is  often  singular  and  hence  non-invertible.  In  order  to  prevent  the 
singularity,  the  piston  of  the  DM,  or  average  actuator  position,  is  set  to  a  fixed  value, 
p.  Equation  (2.12)  can  then  be  written  in  augmented  form  as 


y 

r 

p 

l/N  ■■■  1/N 

(2.14) 


With  this  modification,  the  matrix  is  now  invertible  and  the  phase  of  the  wavefront 
can  be  reconstructed  from  the  slope  measurements.  Other  sensors,  such  as  the  self- 
referencing  interferometer,  which  measure  the  held  directly  instead  of  the  wavefront 
slopes,  do  not  require  this  type  of  phase  reconstruction. 


2. 1.3. 2  Self- Referencing  Interferometer.  Another  type  of  wavefront 
sensor,  known  as  a  self- referencing  interferometer  (SRI  WFS),  directly  measures  the 
wavefront  held  [18].  An  SRI  WFS  works  by  splitting  an  incoming  optical  held  into  a 
reference  held  and  a  beacon  held  as  depicted  in  Figure  2.6.  The  interference  pattern 
resulting  from  the  two  beams  traveling  through  equidistant  paths  can  be  used  to 
measure  the  held  of  the  wavefront.  This  is  done  by  passing  the  reference  held  through 
a  single  mode  hber.  The  hber  spatially  hlters  the  beam  and  provides  a  known  held  at 
the  output.  The  reference  path  also  contains  an  optical  amplifier  and  a  phase  shifter. 
The  phase  shifter  controls  the  phase  of  the  reference  and  the  optical  amplifier  increases 
the  amplitude  of  the  interference  fringes  to  enhance  the  measurement  accuracy. 

If  the  beacon  held  is  defined  as  Ab(x,  y)  exp(i(pb(x,  y)),  and  the  reference  held 
is  defined  as  Ar(x,y),  the  intensity  of  the  n-th  phase  shift  at  the  m-th  subaperture 
center  is 

In(xm,ym)  =  //  [Al(x,y)  +  A2r{x,y)  +  2Ab(x,y)Ar(x,y)  cos((j)b  +  9n)]  dxdy 

"  Am 

(2.15) 
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Figure  2.6:  Diagram  showing  how  an  SRI  WFS  creates  an  interference  pattern  [18]. 


where  9n  is  the  phase  shift  applied  to  the  reference  beam  and  (xm,  ym)  is  the  position 
of  the  subaperture.  The  integral  is  performed  over  the  subaperture  area,  Am.  The 
unknown  beacon  held,  Ub(x,y),  can  be  estimated  from  these  intensity  measurements 
by  using  shifts  of  9n  =  [0,  |,  it,  ^].  This  estimate  is, 

U(xm,ym )  .  {  [A  [Xmi  Vm)  -A  (-^m;  2/m)]  T  i  [A  {Xmi  ym  )  A  i^Xmi  1/m)]  } 

The  amplitude  and  phase  estimates  of  the  beacon  held  can  easily  be  computed  as 


-A-b  {%m  1  Vm) 

—  A  a  V^A^mj  Z/m) 

“1“  [-^4(^777.7  ?/m) 

1 

to 

(2.16) 

7  Urn) 

=  arctan 

-^4(^77 1')  Vm) 

i 

^  h. 

1 

(2.17) 

1\  ym  J 

-*3  Vm  )  _ 

In  this  case,  the  phase  estimate,  0,  lies  between  —7 r  and  7r. 

An  SRI  can  be  designed  to  apply  the  four  phase  shifts  in  a  number  of  ways.  The 
hrst  method  is  to  use  the  hber  phase  shifter  to  apply  each  phase  shift  temporally.  In 
this  case,  four  successive  camera  frames  are  used  to  create  a  single  SRI  measurement. 
Another  method  is  to  divide  the  reference  beam  into  four  separate  beams  and  apply 
the  appropriate  phase  shift  to  each  beam.  The  resulting  interference  intensities  can 
then  be  imaged  simultaneously  using  either  four  individual  detectors  or  four  regions 
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of  a  single  detector.  It  is  also  possible  to  use  a  hybrid  approach,  which  splits  the 
beam  into  two  paths  and  uses  two  successive  frames  from  each  path  to  obtain  the  full 
set  of  four  interference  patterns  necessary  for  a  single  SRI  measurement.  There  are 
advantages  and  disadvantages  to  both  spatial  and  temporal  phase  shifting  techniques. 
When  using  spatial  phase  shifting,  the  same  beacon  field  is  used  for  each  measure¬ 
ment  which  is  advantageous  with  rapidly  evolving  turbulence.  However,  when  using 
multiple  propagation  paths  for  the  different  beams  and  different  detectors  for  each 
interference  pattern,  proper  calibration  and  alignment  become  more  difficult.  Tem¬ 
poral  phase  shifting  avoids  these  calibration  and  alignment  issues  and  also  provides 
a  higher  spatial  resolution  for  a  given  detector,  since  it  is  not  necessary  to  partition 
the  detector  array.  The  drawback  to  temporal  phase  shifting  is  that  measurement 
bandwidth  is  reduced  by  a  factor  of  four  [18]. 

With  either  spatial  or  temporal  phase  shifting,  the  SRI  WFS  provides  a  mea¬ 
surement  of  the  wavefront  phase.  This  measurement  is  used  to  determine  the  shape 
of  the  DM  needed  to  conjugate  the  phase  of  the  wavefront.  The  AO  control  algorithm 
describes  how  the  WFS  measurement  is  used  to  command  the  DM. 

2.1.4  AO  Control.  Although  the  DMs  used  in  AO  have  a  large  number 
of  actuators,  typical  controllers  operate  under  the  assumption  that  each  actuator  is 
independent  and  can  be  controlled  as  a  single-input  single-output  (SISO)  system. 
Typically,  the  same  control  law  is  used  for  each  actuator,  although  some  research  has 
been  done  using  non-uniform  gain  controllers  [21].  The  discrete-time  controller  used 
by  B  archers  is 

c(n  +  1)  =  ac(n)  +  (1  —  a)KDCe(n  +  1)  (2-18) 

where  c{n  +  1)  is  the  next  DM  command,  c(n)  is  the  current  mirror  command,  e(n) 
is  the  residual  phase  measurement  and  a  and  Kdc  are  design  parameters  [3].  The 
parameter  a  is  known  as  the  integrator  gain  and  K Dq  is  the  DC  gain  of  the  controller. 
Equation  (2.18)  is  often  simplified  by  introducing  a  new  parameter,  b,  known  simply 


16 


as  the  gain  which  is  defined  as 


b  —  (1  —  a)Kjjc 


(2.19) 


With  this  substitution,  the  AO  control  law  can  be  written  as 


c(n  +  1)  =  ac(n )  +  be(n  +  1) 


(2.20) 


The  equation  can  also  be  represented  by  its  z  transform 


F(z)  = 


C(z )  bz 


E{z )  z  —  a 


(2.21) 


Figure  2.7  shows  a  diagram  of  this  control  law  used  in  a  typical  AO  system  which  is 
designed  to  regulate  the  residual  phase. 


2.2  Woofer- Tweeter  AO 

One  type  of  AO  system  that  differs  from  the  system  presented  in  Section  2.1, 
is  a  two  mirror  system  known  as  woofer-tweeter  AO.  The  simplest  woofer-tweeter 
arrangement  replaces  the  DM  in  Figure  2.1  with  two  mirrors  in  series.  One  mirror, 
the  woofer,  has  low  bandwidth  and  large  actuator  stroke.  The  other  mirror,  the 
tweeter,  has  high  bandwidth  and  low  actuator  stroke.  The  terminology  is  borrowed 
from  the  familiar  two-speaker  audio  system. 


Figure  2.7:  Typical  AO  controller. 
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2.2.1  Stroke  Requirements.  In  order  to  understand  the  need  for  two  mirrors, 
it  is  necessary  to  estimate  the  actuator  stroke  requirement  of  a  DM  by  finding  the 
expected  phase  variation  over  the  aperture.  The  maximum  phase  variation  can  be 
found  by  computing  the  expected  value  of 


(b(r  +  s)  -<^(r)]2) 


(2.22) 


over  the  aperture,  where  </?(r)  is  the  atmospheric  phase  at  spatial  coordinate  r,  and  s 
is  the  spatial  separation  [10] .  From  this  equation,  an  approximate  stroke  requirement 
can  be  calculated  for  relatively  severe  turbulence.  The  RMS  stroke  requirement  for 
an  aperture  of  75  cm,  r0  of  2  cm,  and  a  wavelength  of  A  is  2.24A  OPD. 

To  avoid  saturation  of  the  DM  actuators,  the  stroke  requirement  should  be  about 
three  times  the  RMS  value,  or  approximately  6.71A  OPD  for  wavelengths  on  the  order 
of  a  micrometer.  A  high-bandwidth  mirror  however,  typically  has  a  stroke  capable 
of  only  ±2  //m  OPD  and  by  itself,  it  is  unable  to  satisfy  the  stroke  requirements  [4], 
Using  two  mirrors  is  one  way  to  achieve  the  necessary  requirement. 


2.2.2  Woofer-Tweeter  Control.  When  controlling  two  mirrors,  the  typical 
control  approach,  described  in  Section  2.1.4  cannot  be  used  directly.  The  additional 
mirror  requires  that  the  controller  correctly  partition  the  mirror  commands  between 
the  two  DMs.  One  approach  to  controlling  a  woofer-tweeter  AO  system  is  to  initially 
present  the  entire  correction  to  the  tweeter  mirror  controller  [4] .  The  tweeter  actuator 
commands  are  then  directed  onto  the  woofer  by  means  of  an  off-loading  matrix  which 
maps  the  tweeter  commands  into  corresponding  woofer  commands.  Because  the  tem¬ 
poral  response  of  the  woofer  is  slower  than  that  of  the  tweeter,  the  high-frequency 
components  of  the  tweeter  commands  are  filtered  by  the  woofer.  The  net  result  of 
the  combined  off-loading  and  woofer  dynamics  is  that  the  woofer  corrects  for  the 
low-frequency  content  while  the  tweeter  corrects  the  high-frequency  content. 
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Figure  2.8:  Woofer-tweeter  AO  controller  with  tweeter  commands  mapped  onto  the 
woofer  [4], 

Another  difference  between  the  woofer  and  tweeter  DM  is  the  spacing  of  the 
actuators.  The  woofer  DM  has  a  lower  actuator  density  than  the  tweeter  mirror. 
Due  to  the  difference  in  the  actuator  density  between  the  two  mirrors,  the  tweeter 
corrects  for  the  high  spatial-frequency  content  while  the  woofer  corrects  for  the  low 
spatial-frequency  content.  Figure  2.8  shows  a  block  diagram  of  this  type  of  control. 
The  offload  matrix,  Ht2Wi  represents  the  mapping  from  the  tweeter  command  to  the 
woofer  commands.  The  addition  of  the  limiters  in  the  block  diagram  is  used  to  prevent 
windup  of  the  integrators  when  the  commanded  correction  exceeds  the  physical  stroke 
limits  of  the  mirror  actuators. 

Other  methods  for  controlling  a  woofer-tweeter  AO  system  use  filters  to  partition 
the  DM  commands  between  the  two  mirrors  [6,7,11].  In  these  controllers,  the  response 
of  the  woofer  mirror  is  included  in  the  tweeter  control  loop  in  order  to  calculate  the 
appropriate  division  of  DM  commands  between  the  woofer  and  the  tweeter  mirrors. 
Thus,  the  tweeter  has  two  control  loops,  one  for  the  low-frequency  content  and  one 
for  the  high-frequency  content. 
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Regardless  of  the  controller  used,  the  purpose  of  the  woofer-tweeter  control  al¬ 
gorithm  is  to  allow  the  dual-mirror  system  to  achieve  the  performance  of  a  system 
which  uses  a  high-bandwidth  mirror  with  no  stroke  limits.  Woofer-tweeter  control 
algorithms  constitute  one  of  the  major  components  of  this  research  effort.  The  other 
major  component,  adaptive  control  of  AO  systems,  is  discussed  in  the  following  sec¬ 
tion. 

2.3  AO  and  Adaptive  Control 

Another  type  of  AO  system  that  differs  from  the  traditional  approach  discussed 
in  Section  2.1  is  an  AO  system  which  uses  an  adaptive  controller.  In  a  traditional, 
or  fixed-gain  controller,  the  control  law  is  programmed  into  the  control  computer  and 
remains  fixed  throughout  the  entire  operation  of  the  AO  system.  In  many  systems, 
particularly  in  laboratory  environments,  an  operator  can  adjust  parameters  in  the 
control  law,  such  as  the  constants  a  and  b  in  Equation  (2.20),  to  achieve  the  desired 
performance  for  a  given  scenario  [3] .  However,  even  when  the  operator  can  change  the 
control  parameters  prior  to  operating  the  AO  system,  they  still  remain  constant  while 
the  system  is  in  use.  In  contrast,  an  AO  system  which  uses  an  adaptive  controller  has 
algorithms  in  place  to  dynamically  vary  the  control  parameters  in  order  to  consistently 
achieve  the  desired  performance,  without  additional  tuning  by  the  operator.  This 
section  discusses  adaptive  controllers  and  their  application  to  AO. 

2.3.1  Adaptive  Control.  While  all  feedback  controllers  react  or  adapt  to 
changes  in  the  environment  through  the  feedback  path,  they  are  not  necessarily  con¬ 
sidered  adaptive.  One  definition  of  an  adaptive  controller  is,  “a  controller  with  ad¬ 
justable  parameters  and  a  mechanism  for  adjusting  the  parameters”  [1].  Because 
adaptive  controller  parameters  vary  in  response  to  the  process  measurements,  they 
are  a  form  of  nonlinear  control.  Adaptive  controllers  are  most  useful  in  cases  where 
the  process  dynamics  or  disturbance  characteristics  change  over  time.  In  such  cases, 
fixed-gain  controllers,  designed  under  the  assumption  that  the  underlying  process  is 
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linear  and  time- invariant,  may  not  perform  adequately.  One  type  of  adaptive  con¬ 
troller  is  a  model  reference  adaptive  system  (MRAS). 

2.3.2  Model  Reference  Adaptive  System.  A  model-reference  adaptive  sys¬ 
tem  is  an  adaptive  control  technique  which  consists  of  a  model  that  describes  the 
desired  response  of  the  regulated  plant,  a  regulator  which  controls  the  plant,  and  an 
adjustment  mechanism  which  varies  the  regulator  parameters.  The  goal  of  the  MRAS 
system  is  to  adjust  the  regulator  parameters  such  that  the  regulated  plant  has  the 
same  behavior  as  the  model.  Figure  2.9  shows  a  block  diagram  of  an  MRAS.  One 
approach  to  designing  an  MRAS  is  the  gradient  method.  In  the  gradient  method,  the 
regulator  parameters  are  adjusted  in  the  direction  that  drives  the  error  between  the 
desired  response  and  the  true  response  to  zero. 

Astrom  and  Wittenmark  outlined  the  steps  of  implementing  an  MRAS  using 
the  gradient  method  [1].  The  first  step  in  the  gradient  method  approach  is  to  define 
a  cost  function,  J,  that  the  controller  tries  to  minimize.  Let 

m  =  y 2  (2.23) 


Figure  2.9:  Model  Reference  Adaptive  System  (MRAS),  is  a  form  of  adaptive  con¬ 
trol  which  finds  the  regulator  parameters  necessary  to  drive  the  plant  such  that  the 
closed-loop  system  behaves  like  the  model  system  [1], 
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where  9  is  the  adjustment  parameter,  and  e  is  the  error  between  the  model  output, 
ym,  and  the  true  system  output,  y.  That  is, 


e  =  y  -ym 


(2.24) 


Note  that  the  adjustment  parameter  9  in  this  derivation  follows  standard  MRAS 
notation  and  is  not  related  to  the  phase  shift  described  in  Section  2. 1.3.2.  With  the 
cost  function  and  error  defined,  the  parameters  are  adjusted  in  the  direction  of  the 
negative  gradient  of  J  in  order  to  minimize  the  cost.  If  the  parameters  of  the  true 
system  are  assumed  to  be  constant,  the  adjustment  of  the  estimated  parameters,  9, 
can  be  written  as 


d,9  dJ  de 

dt  =  ~^~d9  =  ~ie~89 


(2.25) 


where  7  is  an  arbitrary  step  size.  This  adjustment  rule  is  known  as  the  MIT  rule. 
By  using  the  MIT  rule,  the  adaptive  controller  can  determine  the  parameter  values 
necessary  to  drive  the  output  of  the  plant  to  match  the  desired  model  response. 
It  is  important  to  note  that  in  the  MIT  rule  the  adjustment  in  Equation  (2.25)  is 
proportional  to  the  error,  e.  This  means  that  parameter  convergence  slows  once  the 
system’s  performance  begins  to  match  the  model  system.  As  a  result,  the  parameters 
may  not  converge  to  the  true  values.  However,  because  the  goal  of  the  MRAS  is  to 
reduce  e,  it  does  not  matter,  from  a  control  standpoint,  whether  the  parameters  have 
converged.  One  drawback  to  the  MRAS  system  described  here  is  the  fact  that  it  does 
not  guarantee  stability.  A  poor  choice  of  the  parameter  7  can  lead  to  instability  [1]. 
This  limitation  is  one  of  the  reasons  that  lattice  filter  controllers  are  also  be  considered 
in  this  research. 


2.3.3  Lattice  Filter.  Another  type  of  adaptive  control  is  a  lattice-filter-based 
approach  which  has  been  successfully  used  in  several  AO  applications  [8,9,14-17,19]. 
The  controllers  all  use  a  finite  impulse  response  (FIR)  lattice  filter.  FIR  filters  which 
have  no  feedback  consist  of  a  weighted  sum  of  delayed  input  values.  This  structure 
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Figure  2.10:  Typical  structure  of  a  transverse  FIR  filter. 
Stage  1  Stage  M 


Figure  2.11:  A  lattice  filter  consists  of  a  cascade  M  stages  with  an  identical  struc¬ 
ture.  Each  stage  is  characterized  by  its  reflection  coefficient,  nm .  The  name  lattice 
filter  comes  from  the  lattice-like  appearance  of  the  stages. 


guarantees  bounded-input  bounded-output  (BIBO)  stability.  A  lattice  filter  differs 
from  the  more  familiar  transverse  FIR  filters  in  form.  It  can  be  shown  that  any  given 
FIR  transfer  function  which  can  be  implemented  as  a  transverse  filter  can  also  be 
implemented  as  a  lattice  filter  [13].  A  transverse  FIR  filter  has  the  simple  structure 
shown  in  Figure  2.10.  A  lattice  filter  on  the  other  hand,  has  a  more  complex  structure 
as  shown  in  Figure  2.11.  The  name  lattice  filter  refers  to  the  lattice-like  appearance  of 
the  filter  connections.  An  Mth-order  lattice  filter  consists  of  M  stages  with  identical 
structures,  with  each  stage  characterized  by  its  reflection  coefficient,  Krn.  Each  stage 
has  two  inputs  and  two  outputs  that  are  described  by  the  equations 


fm{n)  =  fm- l(n)  +  Kmbm-iin  -  1) 


(2.26) 


and 


bm(n)  =  bm-i(n  -  1)  +  Kmfm-i(n) 


(2.27) 
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Figure  2.12:  A  joint-process  estimator  lattice  filter  consists  of  a  lattice  filter  with 
the  addition  of  weights  which  operate  on  the  backward  prediction  errors. 


with  initial  condition 

fo(n)  =  b0(n )  =  u(n)  (2.28) 

The  quantity  fm(n )  is  known  as  the  m-th  forward  prediction  error  and  represents  the 
difference  between  the  input  u(n)  and  the  predicted  value  based  on  the  m  previous 
inputs  u(n  —  1), . . . ,  u(n  —  m).  Similarly,  the  quantity  bm(n)  is  known  as  the  backward 
prediction  error  and  is  the  difference  between  u{n  —  m)  and  the  m  future  inputs 
u(n), . . .  ,u(n  —  m  +  1). 

By  augmenting  the  structure  shown  in  Figure  2.11  with  an  additional  set  of 
weights,  which  operate  on  the  backward  prediction  errors,  a  lattice  filter  can  be  used 
to  predict  a  correlated  sequence.  This  is  known  as  a  joint-process  estimator  and  is 
shown  in  Figure  2.12.  When  the  filter  input  is  a  correlated  stationary  process,  the 
backward  prediction  errors  have  the  important  property  of  being  orthogonal.  This 
property  of  orthogonal  prediction  errors  allows  the  lattice  filter  joint-process  estimator 
to  adapt  faster  than  a  transversal  filter  because  there  are  no  interactions  between  the 
filter  weights  [13].  In  addition  to  faster  adaptation,  lattice  filters  are  more  numerically 
robust  than  transversal  filters  and  their  modular  structure  allows  the  filter  order  to 
be  increased  without  affecting  previously  calculated  coefficients  [12,13]. 
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Figure  2.13:  Control  architecture  commonly  used  in  AO  controllers  which  use  lattice 
filters  to  provide  additional  input  to  the  traditional  fixed-gain  controller. 

The  previously  mentioned  AO  systems  with  lattice  filter  adaptive  controllers 
have  all  used  the  structure  shown  in  Figure  2.13,  with  the  adaptive  controller  providing 
additional  input  to  the  fixed-gain  controller  [8,9,14-17].  In  the  figure,  the  model  is 
used  to  provide  an  estimate  of  the  residual  phase  measurement,  ym,  based  on  the 
adaptive  controller  output,  ua.  The  signal  ym  is  used  to  calculate  the  error  signal, 
e,  between  the  model  system  and  the  true  system  as  in  Equation  (2.24).  This  error 
signal  is  the  input  to  the  lattice  filter,  L(z).  The  tuning  signal  to  the  lattice  filter  is 
the  true  system  output,  y.  The  gains  of  the  lattice  Elter,  L(z),  are  calculated  such 
that  the  adaptive  control  signal,  ua,  drives  the  residual  phase  measurement,  y ,  to 
zero. 


2.3.4  Modal  Control.  When  dealing  with  a  large  number  of  actuators, 
the  computational  burden  can  become  a  concern,  especially  when  using  complicated 
algorithms  such  as  adaptive  lattice  filters.  One  common  method  of  reducing  the 
computational  burden  of  the  controller  is  to  define  a  set  of  orthogonal  mirror  modes 
and  apply  the  control  law  to  each  individual  mode  instead  of  each  actuator.  As  long 
as  the  number  of  modes  is  less  than  the  number  of  actuators,  then  the  computational 
burden  is  decreased.  Liu  and  Gibson  defined  a  modal  decomposition  that  results  in 
orthogonal  modes  that  are  ordered  according  to  their  spatial-frequency  content  [15]. 
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This  spatial-frequency  ordering  is  useful  because  it  effectively  performs  lowpass  spatial 
filtering  of  the  DM  commands  when  the  higher-order  modes  are  discarded.  This 
can  help  mitigate  high-frequency  noise  and  limit  the  effects  of  mis- registration  [15]. 
The  process  for  determining  these  frequency  weighted  modes  consists  of  two  steps: 
first  hireling  modes  which  are  orthogonal  in  actuator-space  and  then  finding  modes 
which  are  orthogonal  in  both  actuator  and  sensor-space.  In  this  process  the  frequency 
weighting  is  achieved  through  the  use  of  a  lowpass  spatial  filter. 

The  frequency-weighted  modal  decomposition  procedure  requires  only  the  ac¬ 
tuator  geometry  and  a  lowpass  filter.  The  matrices  T  and  S  describe  the  actuator 
geometry  and  the  matrix  F  describes  the  lowpass  filter.  The  columns  of  T  form  a  set 
of  basis  vectors  in  the  space  of  the  master  actuators,  while  S  is  the  slaving  matrix  that 
maps  the  master  actuators  commands  to  commands  for  all  actuators.  The  columns 
of  T  can  be  chosen  such  that  they  represent  DM  shapes  with  no  mean  displacement 
or  tilt.  If  these  constraints  are  not  required,  T  can  be  the  identity  matrix.  Likewise, 
if  there  is  no  slaving,  S  can  be  the  identity  matrix.  F  can  be  the  matrix  representa¬ 
tion  of  any  arbitrary  lowpass  filter.  The  frequency-weighted  modes  are  defined  as  the 
columns  of  the  matrix 

V  =  TX  (2.29) 

where  X  is  the  m  x  m  eigenvector  matrix  which  solves 

(. FST)*(FST)X  =  (ST)*{ST)X  A2  (2.30) 

and  where  A2  is  the  diagonal  matrix  of  the  eigenvalue  magnitudes  squared,  A2,  and  * 
represents  the  conjugate  transpose.  If  only  n  modes  are  needed,  then  only  the  first  n 
columns  of  V  are  used  to  define  the  modes. 

The  modes  described  by  V  are  designed  to  be  orthogonal  in  actuator  space 
and  are  not,  in  general,  orthogonal  in  sensor  space.  In  order  to  find  modes  that  are 
orthogonal  in  both  actuator  and  sensor-space,  the  poke  matrix,  T,  which  describes  the 
relationship  between  actuator  commands  and  sensor  measurements,  is  needed.  The 


26 


new  modes  are  defined  as  the  columns  of  the  matrix 


V  =  TX  (2.31) 

where  X  is  the  n  x  n  eigenvector  matrix  which  solves 

(r  sv)*{vsv)x  =  (sv)*(sv)x  A2  (2.32) 

Using  this  method  of  defining  modes,  which  are  orthogonal  in  both  actuator  and 
sensor  space,  a  lattice  filter  woofer-tweeter  controller  can  now  be  designed.  Building 
on  the  material  covered  in  this  chapter,  Chapter  Ill  discusses  the  design  process  of 
both  MR.AS  and  lattice  filter-based  adaptive  woofer-tweeter  controllers. 
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III.  Simulation  Environment  and  Adaptive  Controller 

Design 


This  chapter  takes  the  topics  in  Chapter  II  and  discusses  the  specikc  details  of 
how  they  can  be  combined  to  create  a  new  adaptive  optics  system  with  improved 
performance  compared  to  existing  approaches.  The  discussion  consists  of  two  parts: 
the  setup  of  the  simulation  and  the  design  of  the  adaptive  controllers.  The  simulation 
details  are  presented  prior  to  developing  the  specific  control  architectures  because 
some  aspects  of  the  design  process,  such  as  defining  the  mirror  modes,  are  dependent 
on  the  type  of  hardware  being  simulated. 

3. 1  Simulation 

In  order  to  design  and  test  a  control  algorithm  for  an  adaptive  woofer-tweeter 
AO  system  it  is  necessary  to  have  a  working  simulation  of  all  of  the  major  components 
of  an  AO  system.  As  discussed  in  Section  2.1,  these  components  are  the  atmospheric 
turbulence,  the  DM,  the  WFS,  and  the  control  algorithm.  Several  commercial  software 
packages  exists  for  simulating  these  components.  Because  they  were  readily  available 
and  relatively  easy  to  learn,  the  MATLAB  packages  WaveProp  and  AOTools  were 
used  for  these  simulations  [5]. 

3.1.1  Atmosphere.  One  way  of  generating  atmospheric  turbulence  with  the 
desired  statistics  is  to  apply  a  phase  screen  to  a  simulated  complex  optical  field.  The 
phase  screen  adds  or  subtracts  phase  from  the  optical  field  to  simulate  the  effects  of 
atmosphere  on  a  wavefront.  Figure  3.1  shows  an  example  of  a  phase  screen  realization 
using  the  Kolmogorov  turbulence  model  given  in  Equation  (2.3),  with  r0  =  2  cm.  By 
applying  a  phase  screen  with  the  appropriate  value  of  r0  to  the  complex  optical  field, 
it  is  possible  to  simulate  the  desired  turbulence  to  test  the  adaptive  AO  controller. 
While  a  phase  screen  provides  the  correct  spatial  characteristics  of  the  turbulence, 
it  is  also  necessary  to  simulate  the  temporal  variations  of  the  turbulence.  Given  a 
wind  velocity,  V,  the  corresponding  Greenwood  frequency,  fa,  can  be  determined 


Kolmogorov  Phase  Screen  for  rQ  =  2  cm 
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Figure  3.1:  An  example  of  a  phase  screen  generated  using  the  Kolmogorov  turbu¬ 
lence  model  with  r0  =  2  cm. 

from  Equation  (2.11).  To  simulate  the  temporal  effects,  the  phase  screen  is  simply 
shifted  laterally  by  a  distance  V  At,  where  At  is  the  simulation  time  step. 

Because  the  focus  of  this  research  is  to  develop  an  adaptive  control  algorithm, 
only  a  single  phase  screen  is  used  and  all  the  EM  fields  are  simulated  in  the  pupil 
plane  of  the  telescope.  This  approach  eliminates  the  need  for  additional  elements  in 
the  simulation,  such  as  telescopes  which  crop  and  demagnify  the  held,  and  computa¬ 
tionally  expensive  wave  propagations  which  simulate  both  scintillation  of  the  beacon 
wavefront  and  the  directed  energy  performance  at  the  target  location.  Neither  the 
telescope  nor  the  full-path  wave  propagations  are  necessary  for  demonstrating  the 
feasibility  of  a  new  AO  control  architecture.  The  inclusion  of  either  the  telescope  or 
propagations  would  not  change  the  design  nor  implementation  of  the  adaptive  control 
algorithms. 

3.1.2  Simulation  Parameters.  When  simulating  the  performance  of  the 
adaptive  woofer-tweeter  AO  system,  two  different  turbulent  scenarios  are  considered. 
The  two  scenarios  are  the  same  ones  used  by  Brennan  and  Rhoadarmer  [4].  The  first 
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Table  3.1:  Simulation  Parameters 


Parameter 

Low  Turbulence 

High  Turbulence 

r0 

6  cm 

2  cm 

V 

20  m/s 

20  m/s 

/g 

142  Hz 

427  Hz 

scenario  consists  of  fairly  correlated  turbulence,  which  should  be  easily  compensated 
by  the  standard  fixed-gain  AO  system.  The  second  scenario  consists  of  less-correlated 
turbulence  which  is  more  challenging  for  the  standard  controller.  It  is  expected  that 
the  adaptive  controller  can  improve  the  performance  of  the  AO  system  in  more  dif¬ 
ficult  turbulence  while  maintaining  the  performance  of  the  AO  system  in  the  case  of 
the  lower  turbulence.  Table  3.1  lists  the  simulation  parameters  of  the  two  different 
scenarios. 

3.1.3  Deformable  Mirrors.  As  with  the  turbulent  scenarios,  the  woofer  and 
tweeter  DM  specifications  are  taken  from  Brennan  and  Rhoadarmer  [4].  Both  the 
woofer  and  the  tweeter  DMs  are  set  equal  to  the  simulated  aperture  size  of  76  cm. 
The  tweeter  DM  has  39  actuators  across  the  aperture  with  a  spacing  of  2  cm  for  a 
total  of  1521  actuators  over  the  square  aperture.  The  woofer  DM  has  20  actuators 
across  the  aperture  with  a  spacing  of  4  cm  for  a  total  of  400  actuators.  To  simulate 
a  continuous  facesheet  DM,  the  surface  of  each  mirrors  is  approximated  by  a  bilinear 
interpolation  of  its  actuator  commands.  As  shown  in  Figure  3.2,  the  two  mirrors  are 
positioned  such  that  each  woofer  actuator  is  aligned  with  a  corresponding  tweeter 
actuator. 

In  order  to  simulate  a  circular  aperture,  only  some  of  the  actuators  over  the 
square  DM  are  active.  Figure  3.3  shows  the  layout  of  the  active  actuators  for  both 
the  woofer  and  the  tweeter  DMs.  The  total  number  of  active  actuators  is  1245  for  the 
tweeter  and  316  for  the  woofer.  The  alignment  shown  in  Figure  3.2  makes  it  easy  to 
determine  the  offloading  matrix  Ht2w ■  If  the  active  actuators  in  the  39  x  39  tweeter 
grid  are  represented  by  a  1245  x  1  vector  and  the  woofer  is  similarly  represented  by 
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Woofer  Actuator 
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Figure  3.2:  Woofer-tweeter  actuator  alignment.  Each  woofer  actuator  corresponds 
to  a  tweeter  actuator.  The  spacing  between  woofer  actuators  is  twice  that  of  the 
tweeter  actuators  giving  the  tweeter  a  higher  spatial  resolution.  Not  all  of  the  actua¬ 
tors  are  depicted  here. 

a  316  x  1  vector,  Ht2W  is  a  316  x  1245  matrix  which  maps  the  tweeter  actuators  to 
the  corresponding  woofer  actuators.  Each  tweeter  actuator  that  has  a  corresponding 
woofer  actuator  is  represented  by  a  one  in  the  matrix.  All  other  entries  in  the  matrix 
are  zeros.  In  other  words,  Ht2w  effectively  down-samples  or  decimates  a  vector  of 
tweeter  commands  to  produce  the  corresponding  set  of  woofer  commands.  The  other 
matrix,  Hw2t,  is  used  to  determine  the  effect  of  the  woofer  DM  at  the  tweeter  DM 
actuator  locations,  in  order  to  accurately  combine  the  correction  of  both  mirrors. 
This  matrix  is  not  needed  explicitly  since  it  is  automatically  accounted  for  in  the 
simulation  when  the  mirror  surface  is  applied  to  the  optical  path.  If  the  matrix  was 
needed,  it  would  simply  be  the  matrix  representation  of  a  bilinear  interpolation  of  the 
woofer  DM  commands  to  the  higher  resolution  of  the  tweeter  DM  commands. 

In  the  simulations,  the  tweeter  DM  is  assumed  to  have  an  infinitely  fast  response 
time,  compared  to  the  sensor  rate,  and  can  be  modeled  by  the  transfer  function 

T(z)  =  1  (3.1) 


31 


0.3 


Figure  3.3:  To  simulate  a  circular  aperture  only  some  of  the 
actuators  are  active.  The  total  number  of  active  tweeter  actua¬ 
tors  is  1245  as  shown  in  (a).  The  total  number  of  active  woofer 
actuators  is  316  as  shown  in  (b). 


Table  3.2:  Woofer  and  Tweeter  DM  Parameters 


Parameter 

Tweeter  DM 

Woofer  DM 

Integrator  gain 

dr  —  0.99 

aw  =  0.99 

Gain 

bx  —  0.4 

b\y  =  0.1 

Stroke  Limit 

Lt  =  ±2.5/im 

Lw  =  ±6.5/im 

Subaperture  Width 

d,T  =  2  cm 

dw  =  4  cm 

On  the  other  hand,  the  woofer  has  a  slower  response  which  can  be  modeled  by  the 
transfer  function 

W(z)  =  (3.2) 

The  tweeter  mirror  has  a  maximum  stroke  limit  of  ±2.5/im  OPD  while  the  woofer  DM 
has  a  maximum  stroke  limit  of  ±6.5/im  OPD.  Each  mirror  is  driven  by  the  standard 
controller  given  by  Equation  (2.20).  For  the  given  mirrors,  the  optimal  values  for 
the  control  law  were  determined  by  Brennan  and  Rhoadarmer  [4].  The  values  for  the 
control  laws,  as  well  as  a  summary  of  both  DM  characteristics  are  given  in  Table  3.2. 
The  woofer  and  tweeter  DM  gains  differ  due  to  the  dissimilar  temporal  responses  of 
the  mirrors. 
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3.1.4  Wavefro7it  Sensor.  The  wavefront  sensor  used  in  the  simulations  is 
an  SRI  WFS  with  spatially  applied  phase  shifts.  Each  SRI  WFS  subaperture  corre¬ 
sponds  to  one  tweeter  actuator  giving  a  one-to-one  correspondence  between  the  SRI 
WFS  measurements  and  the  tweeter  DM  commands.  Because  the  SRI  WFS  provides 
wrapped  phase  between  —  tt  and  tt,  a  phase  unwrapping  algorithm  is  needed.  In  the 
simulation  software,  various  phase  unwrapping  algorithms  are  available.  The  common 
least-squares  phase  unwrapping  algorithm  is  inadequate  in  the  high  turbulence  case 
because  the  incidence  of  27r-phase  differences  in  adjacent  subapertures  is  too  great. 
Instead,  the  proprietary  WaveProp  unwrapping  function  xphase  was  used  to  avoid 
these  problems. 

3.2  Adaptive  Controller  Design 

The  basic  approach  to  designing  an  adaptive  controller  for  a  woofer-tweeter 
AO  system  employed  in  this  research  is  to  adaptively  control  the  tweeter  mirror 
only  and  then  transfer  the  average  tweeter  commands  onto  the  woofer  as  described 
in  Section  2.2.2.  This  approach  simplifies  the  design  of  the  controller,  especially 
when  determining  the  modes  for  the  adaptive  lattice  filter.  With  this  approach,  only 
the  tweeter  modes  need  to  be  defined,  and  thus  it  is  not  required  to  find  modes 
that  accurately  describe  two  mirrors  with  different  characteristics.  In  this  section, 
adaptive  control  of  the  tweeter  DM  which  is  then  transferred  onto  the  woofer  DM  is 
implemented  with  both  an  MRAS  controller,  as  shown  in  Section  3.2.1,  and  with  a 
lattice  filter-based  controller,  as  shown  in  Section  3.2.2. 

3.2.1  MRAS  Co?itroller  Design.  Because  an  MRAS  controller  is  relatively 
simple  to  both  implement  and  debug,  it  is  used  as  a  preliminary  design  to  verify  the 
feasibility  of  the  general  adaptive  AO  approach.  The  successful  MRAS  approach  is 
then  implemented  using  the  more  complicated  lattice-filter-based  controller. 

One  problem  that  occurs  when  attempting  to  use  an  MRAS  with  an  AO  system 
is  that  in  AO,  the  commanded  input,  ucmd,  is  always  zero.  As  is  shown,  part  of  the 
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Figure  3.4:  MRAS-based  adaptive  AO  control. 


parameter  adjustment  is  proportional  to  ucmd,  thus,  with  no  commanded  input,  the 
adjustment  cannot  occur.  One  way  to  address  this  problem  is  to  use  the  approach 
employed  by  Gibson  et  al.  [9,14-16]  shown  in  Figure  2.13.  In  this  approach  the  adap¬ 
tive  controller  is  used  in  addition  to  the  standard  fixed-gain  controller.  The  adaptive 
controller  is  given  the  same  input  as  the  fixed-gain  controller,  which  is  the  error  sig¬ 
nal  ucrnd  —  y.  The  controller  then  creates  an  additional  adaptive  control  signal,  which 
appears  as  an  additional  error  to  the  fixed-gain  controller  in  order  to  drive  the  overall 
system  toward  the  desired  response.  Although  Gibson  et  al.  used  a  lattice-filter-based 
adaptive  controller  instead  of  an  MR  AS,  the  general  approach  is  still  the  same.  Fig¬ 
ure  3.4  shows  the  closed-loop  AO  system  with  an  MRAS  adaptive  controller  being 
used  in  this  manner.  When  implemented  in  this  way,  the  model  system  in  the  MRAS 
controller  is  an  approximation  of  the  closed-loop  transfer  function,  Y(z)/Ua(z ),  of  the 
AO  system  using  only  the  traditional  fixed-gain  controller.  Because  the  controller 
is  deterministic  and  known,  the  MRAS,  in  effect,  attempts  to  compensate  for  any 
modeling  errors  or  non-linearities  in  the  DM  and  WFS  as  well  as  changes  in  the 
atmospheric  turbulence. 

In  order  to  implement  the  MIT  rule  as  discussed  in  Section  2.3.2,  it  is  necessary 
to  find  the  sensitivity  derivative,  de/d6,  from  Equation  (2.25).  The  first  step  in 
determining  the  sensitivity  derivative  is  to  define  the  desired,  or  model,  response. 
Because  the  desired  response  of  the  woofer-tweeter  arrangement  is  to  behave  like  an 
ideal  tweeter  with  no  stroke  limits,  the  woofer  is  neglected  when  deriving  the  model 
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system  for  the  MRAS  controller.  The  tweeter  controller  is 


Ct(z) 


bx 

z  —  ax 


and  the  plant  consisting  of  the  DM  and  WFS  is 


(3.3) 


P(z)  =  1 


(3.4) 


Assuming  a  unit  delay  on  the  feedback  path,  the  MRAS  model  transfer  function  is 


— Ct(z)P(z)  —bxz 

1  +  \CT{z)P(z)  z2  -  aTz  +  bT 


(3.5) 


The  next  step  is  to  find  the  true  system  response.  Because  the  true  system,  G(z), 
is  unknown,  it  is  assumed  to  have  the  same  form  as  Equation  (3.5)  with  unknown 
constants  hi,  k2,  and  k 3.  That  is, 


G(z) 


Y{z) 

Ua(z) 


—k\z 

z2  -  k2z  +  fa 


Next,  let  the  output  of  the  MRAS  regulator,  ua  be  defined  as 


(3.6) 


ua{t)  =  t0uc(t)  -  s0y(t  -  1) 


(3.7) 


where  to  and  So  are  the  adjustment  parameters.  Substituting  Equation  (3.7)  into 
Equation  (3.6),  the  transfer  function,  H(z),  from  the  commanded  input  to  the  true 
plant  output  is  found  to  be 


H(z) 


Y» 

(3,8) 

Y(z) 

— kitoz 

(3.9) 

Uc(z ) 

z2  -  k2z  +  (h3  -  s0ki) 
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Finally,  the  sensitivity  derivatives,  de/dto  and  de/ds^  can  be  determined  from 
Equation  (3.9).  Because  the  model  response,  yrn ,  does  not  depend  on  the  adjustment 
parameters,  to  and  so,  it  can  be  ignored  when  taking  the  partial  derivatives.  The 
sensitivity  derivatives  are 


de 
d  t0 
de 
d  s0 


— k\z 


z2  -  k2z  +  (k:i  -  s0ki) 
—t0klz 

(z2  -  k2z  +  (k3  -  s0ki)) 


Uc(z ) 

~^Uc(z) 


(3.10) 

(3.11) 


These  equations  illustrate  the  previously  mentioned  fact  that  without  an  input  signal, 
uc ,  no  adaptation  can  occur  since  the  sensitivity  derivatives  are  identically  zero.  The 
second  equation,  (3.11),  can  be  simplified  and  rewritten  as  a  function  of  Y(z)  using 
the  relationship  given  in  Equation  (3.9).  The  simplified  sensitivity  derivative  is 


de 

ds0 


k\ 


z2  -  k2z  +  (k3  -  s0ki) 


Y(z) 


(3.12) 


Because  the  parameters  of  the  true  system,  £q,  k2 ,  and  k:i,  are  unknown,  the 
adjustment  rule  cannot  be  implemented  using  Equations  (3.10)  and  (3.12).  In  order 
to  implement  the  controller  it  is  necessary  to  make  some  approximations.  The  value 
k\  appears  in  the  numerator  of  both  derivatives.  Because  the  sign  of  k\  is  known, 
it  can  be  incorporated  into  the  arbitrary  parameter  7  by  letting  7  =  7 k\.  As  the 
parameters  converge  on  their  desired  values,  the  denominator  approaches  that  of  the 
desired  response,  z2  —  cltz  +  br-  Thus,  the  adjustment  rule  can  be  based  on  these 
desired  parameters.  With  these  approximations  and  by  applying  Equation  (2.25),  the 
parameter  update  becomes 


dto 

dt 


ye 


z2  —  cltz  +  br 


dso 

dt 


—ye 


1 

—  CLtZ  + 


(3.13) 

(3.14) 
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In  order  to  implement  the  MRAS  controller  with  the  woofer-tweeter  AO  sys¬ 
tem,  the  adjustment  rule  was  applied  to  each  tweeter  actuator  independently.  In  the 
discrete-time  simulation,  the  sensitivity  derivative  terms  were  propagated  at  each  time 
step  using  the  time  domain  representation  of  their  transfer  function.  Both  updated 
parameters,  for  each  actuator,  were  initialized  to  one.  The  parameters  were  updated 
by  simple  Euler  integration  of  Equations  (3.13)  and  (3.14),  that  is 

to(n)  =  tQ(n  -  1)  + 
s0(n)  =  s0(n  -  1)  + 

where  At  is  the  simulation  time  step. 

The  adaptation  parameter,  7,  was  set  to  1015.  While  this  number  seems  very 
large,  it  is  necessary  due  to  the  fact  that  the  sampling  time  is  10-4  seconds,  the  error 
signal,  e,  is  on  the  order  of  ICE8  meters,  and  the  command  signal,  uc.  is  on  the  order 
of  10''.  If  7  is  too  small,  the  adaptation  occurs  too  slowly  to  be  of  any  benefit  [1]. 
With  an  appropriate  step  size,  the  MRAS  controller  adjusts  the  controller  parameters 
and  drives  the  AO  system  to  have  the  desired  performance. 

3.2.2  Lattice  Filter  Controller  Design.  The  lattice-filter-based  adaptive  con¬ 
troller  uses  the  same  general  control  architecture  as  the  MRAS-based  controller,  with 
the  adaptive  controller  providing  additional  input  to  the  fixed-gain  controller.  Specif¬ 
ically,  it  used  the  same  structure  shown  in  Figure  2.13  and  discussed  in  Section  2.3.3. 
Like  the  MRAS-based  controller,  the  lattice  filter-based  controller  also  uses  the  model 
transfer  function  from  Equation  (3.5).  The  lattice  filter  implemented  is  the  gradient 
adaptive  lattice  (GAL)  filter  described  by  Haykin.  The  relevant  equations  are  sum¬ 
marized  in  Equations  (3.17)-(3.27)  [12].  A  more  thorough  treatment  of  the  GAL  can 
be  found  Chapter  12  of  Haykin’s  book  with  prefatory  material  in  Chapters  1,  3,  5,  and 
6  [12].  The  goal  of  the  GAL  is  to  identify  the  reflection  coefficients,  Km,  to  minimize 


(3.15) 

(3.16) 
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the  cost  function 


Jfb,m  ~  i E[\fm(n)\ 2  +  \bm(n)\2] 


(3.17) 


where  fm(n)  and  bm{n )  are  the  forward  and  backward  prediction  errors  mentioned  in 
Section  2.3.3  and  E  is  the  expectation  operator.  Using  Equations  (2.26)  and  (2.27), 
the  derivative  of  Jfb,m  with  respect  to  Km  can  be  found  to  be 


r)  T 

(S[|/™-i(«)|2]  +  E[\bm-i(n  -  1)|2]) 


O  Z7TJ. 


Setting  the  gradient  equal  to  zero  and  solving  for  Km  gives  the  optimal  reflection 
coefficient 


2E\bm-i(n-l)fm-1(n)] 

U[|/m_1(n)|2  +  |6m_1(n-l)|2]  1  J 

Assuming  the  process  is  ergodic,  the  expectation  operator  can  be  replaced  by  time 
averages  to  give  the  time- varying  estimate  of  the  reflection  coefficient,  km{n).  This 
estimate,  based  on  time  averages,  can  be  written  recursively  as 


„  ,  v  „  ,  fm-i(n)bm(n)  +  bm-i(n  ~  l)fm(n) 

Km{n)  =  Km(n  -  1) - - - — - 


(3.20) 


where  Sm-i(n)  is  the  total  energy  of  both  prediction  errors  at  time  n  at  the  ontpnt 
of  stage  m  —  1.  The  total  energy  is  given  by 


£m-i(n)  =  £m-i(n  -  1)  +  |/m_i(n)|2  +  |6m_i(n  -  1)|2  (3.21) 


The  initial  values  for  Equations  (3.20)  and  (3.21)  are  fcm(0)  =  0  and  £m_i(0)  = 
0.  Equation  (3.20)  can  be  rewritten  to  incorporate  a  step-size  parameter,  jl,  which 
leads  to  faster  convergence  when  the  forward  and  backward  prediction  errors,  and 
consequently  £m-u  are  small.  The  modified  equation  is 


km{n) 


„  ,  ^  At  fm-i(n)bm(n)  +  bm_i(n 

Km{n  -  1)  -  - - — - - - - - 

£m—  liP')  £m—  IkE) 


•  1)/, 


myl) 


(3.22) 
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Similarly,  Equation  (3.21)  can  be  modified  with  the  addition  of  the  parameter  f3  to 
allow  for  non-stationary  environments.  The  modified  equation  is 

Sm-i(n)  =  f3£m-i(n  -  1)  +  (1  -  (3)  (|/m_i(n)fs  +  |6m-i (n  -  1)|2)  (3.23) 

where  f3  is  a  scalar  in  the  interval  0  <  /?  <  1  and  serves  to  provide  a  memory  of  the 
past  value  of  £m-i  when  computing  km{n)  [12]. 

The  next  part  of  the  GAL  algorithm  described  by  Haykin  is  the  desired  response 
or  joint-process  estimator,  which  predicts  the  desired  response  dm(n )  based  on  the 
backward  prediction  errors,  that  is 

m 

dm{n)  =  ^2  hk(n)bk(n) 

k= 0 

=  dm_i(n)  +  hm(n)bm(n)  (3.24) 

The  values  of  hm(n )  are  determined  using  a  similar  gradient  approach  that  yields  the 
recursive  time  update  equation 

h'miP'  T  1 )  =  hm(n)  T  tt—  bm(ji)6.m(ji)  (3.25) 

||bm(n)||- 

where  em{n )  is  the  estimation  error  given  by 

em(n)  =  d(n )  -  dm(n)  (3.26) 

The  quantity  b m(n)  is  the  input  vector  of  the  backward  prediction  errors  bo(n),  bi(n), . . . , 
bm(n),  and  fi  is  the  step-size  parameter.  The  Euclidean  norm  of  b m(n)  can  be  recur¬ 
sively  determined  by 

m 

1 1 bjn (n) 1 1 2  =  J2  Mn)\2 

k= 0 

=  |  |bm_i(n)|p  +  \bm(n)\-  (3.27) 
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Once  the  lattice  filter  has  been  designed  according  to  Equations  (3.20)-(3.27), 
the  next  step  is  to  determine  a  set  of  orthogonal  mirror  modes.  Having  the  lattice 
filter  operate  on  a  set  of  orthogonal  mirror  modes  serves  to  both  reduce  number  of 
input  channels  to  the  filters  and  allows  the  filter  to  adapt  to  both  spatial  and  temporal 
variations.  When  the  lattice  filter  is  applied  to  each  actuator  independently,  it  adapts 
to  only  the  temporal  variations  at  each  actuator  location.  However,  when  the  lattice 
filter  is  applied  to  the  orthogonal  mirror  modes,  the  filter  adapts  to  the  temporal 
variations  of  each  mode.  Because  each  mode  contains  spatial  information  from  the 
entire  mirror,  as  the  filter  adapts  to  each  mode  it  is  also  adapting  to  spatial  variations 
across  all  actuators. 

The  modes  were  determined  using  the  procedure  discussed  in  Section  2.3.3.  In 
order  to  determine  the  modes,  it  is  necessary  to  define  the  matrices  F,  S,  and  T  in 
Equation  (2.30).  In  these  simulations,  with  no  slaving  and  no  constraints  on  T,  both 
matrices  S  and  T  are  simply  the  identity  matrix.  The  matrix  F  was  created  using  the 
lowpass  filter  shown  in  Figure  3.5.  The  lowpass  filter  is  the  same  as  the  filter  used  by 
Liu  and  Gibson  [15].  The  filter  was  generated  by  creating  a  one-dimensional  eighth- 
order  FIR  filter,  iq,  with  a  normalized  cutoff  frequency  of  0.5  using  the  MATLAB 
command  hrl.  A  two-dimensional  filter,  F2,  was  created  from  the  one- dimensional 
filter  using  the  MATLAB  command  ftrans2,  which  performs  a  McClellan  transfor¬ 
mation.  The  matrix  F  is  created  such  that  it  maps  a  vector  of  DM  commands  to  a 
vector  of  commands  that  have  been  spatially  filtered  by  F2. 

Once  the  matrices  F,  S,  and  T  have  been  defined,  the  modes  can  be  calculated 
directly  from  Equation  (2.30).  When  applying  Equation  (2.30),  it  is  necessary  to 
sort  the  resulting  eigenvalues  in  descending  order  to  obtain  the  correct  frequency 
ordering  prior  to  selecting  the  n  columns  of  the  matrix  V  which  are  used  in  the  control 
algorithm.  Naturally,  the  columns  of  the  eigenvector  matrix  also  need  to  be  sorted 
to  maintain  their  relationship  with  the  sorted  eigenvalues.  Because  the  sensor  used  is 
an  SRI  WFS,  where  each  phase  measurement  corresponds  directly  with  one  actuator, 
the  actuator-space  and  the  sensor-space  are  identical.  This  eliminates  the  need  for 
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Figure  3.5:  Spatial  lowpass  filter  used  to  determine  the  DM 
modes  used  in  the  lattice-filer  adaptive  controller. 

applying  Equation  (2.32),  due  to  the  fact  that  the  modes  are  already  orthogonal  in 
both  actuator-space  and  sensor-space.  Figure  3.6  depicts  the  first  24  modes  used  for 
these  simulations. 

Now  that  the  modes  have  been  defined,  it  is  necessary  to  discuss  how  the  modes 
are  incorporated  into  the  woofer-tweeter  control  architecture.  Because  the  modes  were 
defined  for  the  tweeter  mirror,  they  can  be  applied  only  to  the  tweeter  portion  of  the 
control  loop.  Figure  3.7  shows  how  this  is  accomplished.  At  the  input  to  the  adap¬ 
tive  controller,  the  SRI  WFS  measurements  are  multiplied  by  V+  to  map  the  phase 
measurements  from  actuator  space  to  modal-space,  where  +  represents  the  pseudoin¬ 
verse.  Before  the  adaptive  controller  output  is  applied  to  the  traditional  controller,  it 
is  mapped  back  into  actuator-space  by  multiplying  by  V.  Thus,  the  tweeter  is  con¬ 
trolled  adaptively  in  modal-space  while  the  traditional  woofer  and  tweeter  controllers 
operate  in  actuator-space  as  shown  in  Figure  3.7. 

With  modes  defined  and  a  control  architecture  which  uses  them  in  a  woofer- 
tweeter  AO  system,  the  design  of  the  adaptive  woofer-tweeter  controller  is  complete. 
Chapter  IV  examines  the  performance  of  the  designed  controllers  using  the  simulation 
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First  24  DM  modes  used  ii 

i  the  lattice  filter  adaptive  con- 

troller.  The  modes  are  orthogonal  and  ordered  according  to  their  spatial- 
frequency  content. 
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parameters  mentioned  in  this  chapter.  The  controller  performance  in  the  simulations 
is  used  to  determine  the  appropriate  controller  settings,  such  as  the  lattice  filter  order 
and  required  training  time,  and  to  evaluate  the  design. 
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IV.  Results  and  Data  Analysis 


This  chapter  presents  the  simulations  of  the  controllers  developed  in  Chapter  III 
and  compares  them  to  a  simulation  of  the  fixed-gain  woofer-tweeter  system. 
The  fixed-gain  woofer-tweeter  system  is  simulated  first  so  that  a  baseline  can  be 
established.  The  other  simulations  are  compared  both  to  the  baseline  and  to  each 
other.  The  second  simulation  uses  the  comparatively  simple  MR.AS  controller  in 
order  to  validate  the  control  architecture  of  adaptively  controlling  only  the  tweeter 
commands  and  then  transferring  them  onto  the  woofer.  Lastly,  the  lattice-filter-based 
controller  is  simulated.  This  lattice  filter  simulation  is  broken  into  two  parts:  one  in 
actuator-space  and  one  in  modal-space.  The  simulation  in  actuator-space  is  used  to 
determine  the  required  training  time  of  the  filter  and  the  necessary  filter  order.  The 
modal-space  simulation  is  then  performed  using  the  same  training  time  and  filter 
order.  Thus,  simulations  are  presented  in  order  of  increasing  complexity  and  are  used 
as  stepping  stones  to  the  next  level  of  complexity.  This  approach  mirrors  the  design 
process  that  led  to  the  development  of  the  lattice  filter  controller  operating  on  the 
system  modes. 

4-1  Woofer- Tweeter  Simulation  Results 

In  order  to  assess  the  performance  of  the  adaptive  controllers  presented  in  Chap¬ 
ter  III,  it  is  necessary  to  establish  a  baseline.  This  was  done  by  simulating  a  woofer- 
tweeter  AO  system  that  uses  traditional  fixed-gain  controllers.  This  is  the  same  system 
which  is  augmented  with  the  adaptive  control  loop  in  subsequent  simulations.  The 
first  step  in  ensuring  that  the  woofer-tweeter  simulation  is  working  is  to  verify  that 
the  atmospheric  correction  is  being  properly  partitioned  between  the  two  mirrors. 
Figures  4.1  through  4.3  form  a  time  sequence  showing  the  corrected  turbulence,  the 
woofer  and  tweeter  commands  used  to  achieve  the  correction,  as  well  as  the  residual 
phase  measurement  from  the  SRI  WFS.  As  seen  in  Figure  4.1  the  initial  correction  is 
presented  entirely  to  the  tweeter  mirror.  The  woofer  initially  has  no  command  due 
to  the  latency  of  the  transferring  procedure.  Figure  4.2  shows  that  after  0.5  ms  the 
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tweeter  commands  are  beginning  to  be  transferred  onto  the  woofer.  However,  the 
tweeter  is  still  supplying  the  majority  of  the  correction  as  seen  by  the  magnitude  of 
the  commands.  Figure  4.3  shows  that  after  2.0  ms  the  woofer  mirror  has  assumed  the 
majority  of  the  low-spatial-frequency  corrections  and  the  tweeter  mirror  is  correcting 
for  the  higher-spatial- frequency  content. 

Another  way  to  visualize  the  transferring  of  the  tweeter  commands  onto  the 
woofer  is  to  pick  a  point  in  space  where  both  mirrors  have  an  actuator  and  look  at 
the  time  history  of  both  actuators  at  that  point.  If  the  woofer-tweeter  controller  is 
working  properly,  the  initial  correction  should  be  handled  exclusively  by  the  tweeter 
actuator.  Over  time  the  woofer  actuator  should  assume  the  large  amplitude  low- 
frequency  component  while  the  tweeter  actuator  corrects  only  the  remaining  smaller 
amplitude  high-frequency  content.  Figure  4.4  shows  that  this  is  indeed  the  case. 
Furthermore,  this  figure  highlights  the  need  for  woofer-tweeter  AO  since  the  total 
correction,  shown  by  the  dashed  black  line,  frequently  exceeds  the  tweeter  actuator’s 
physical  limits  of  ±2. 5/un. 

Figure  4.5  shows  the  performance  of  the  woofer-tweeter  AO  system  as  measured 
by  the  Strehl  ratio  for  both  the  low  turbulence  and  the  high  turbulence.  The  average 
Strehl  ratio  in  the  low  turbulence  case  was  0.9063.  With  the  stronger  turbulence, 
the  average  Strehl  ratio  drops  to  0.4625.  For  comparison  purposes,  a  single  mirror 
case  is  shown  in  which  only  the  tweeter  was  simulated.  As  expected,  the  woofer- 
tweeter  AO  system  outperforms  the  single  mirror  system  in  both  cases.  With  a  single 
mirror,  the  average  Strehl  ratios  are  0.8921  and  0.3725  for  the  low  and  high  turbulence 
cases,  respectively.  The  large  degradation  in  the  single  mirror  performance  in  strong 
turbulence  is  primarily  due  to  saturation  of  the  tweeter  mirror  actuators. 

4-2  MR  AS  Simulation  Results 

With  the  baseline  established,  the  performance  of  the  different  controllers  can 
now  be  assessed.  When  the  MRAS  controller  is  implemented  as  shown  in  Figure  3.4, 
the  injection  of  the  adaptive  control  signal,  ua,  effectively  changes  the  gain  of  the 
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Time  =  0.1  (ms) 


Figure  4.1:  The  upper-left  figure  shows  the  corrected  phase 
of  the  simulated  optical  held,  while  the  upper  right  shows  the 
corresponding  SRI  measurement.  The  lower  right  shows  the 
woofer  DM  command  and  the  lower  left  shows  that  tweeter  DM 
command.  Initially  the  correction  is  handled  exclusively  by  the 
tweeter  due  to  the  latency  in  the  transfering  calculation. 


Time  =  0.5  (ms) 


Figure  4.2:  After  0.5  ms  the  tweeter  commands  are  beginning 
to  be  transferred  onto  the  woofer.  The  tweeter  is  still  providing 
the  majority  of  the  correction. 
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Corrected  Phase 
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Time  =  2.0  (ms) 


Figure  4.3:  After  2.0  ms  the  woofer  is  correcting  most  of  the 
low  spatial-frequency  components.  The  remaining  high  spatial- 
frequency  components  are  corrected  by  the  tweeter. 


Woofer/Tweeter  Demand 


Figure  4.4:  This  plot  shows  the  history  of  a  woofer  actu¬ 

ator  and  a  tweeter  actuator  that  are  positioned  at  the  same 
location  in  space.  The  tweeter  assumes  the  low-amplitude, 
high-frequency  corrections  and  the  woofer  corrects  the  large- 
amplitude,  high-frequency  content. 
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Strehl  Ratio  (r0  =  6  cm) 


Time  (ms) 


(a) 

Strehl  Ratio  (r0  =  2  cm) 


(b) 


Figure  4.5:  Woofer-tweeter  baseline  performance  in  both 

(a)  low  turbulence  and  (b)  high  turbulence.  The  red  line  shows 
the  performance  of  the  woofer-tweeter  AO  system  and  the  blue 
line  shows  the  performance  of  the  single  mirror  system  using  the 
tweeter  specifications. 
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Parameter  so 


Figure  4.6:  Spatial  variation  of  adjustment  parameter  s0. 

system.  Because  the  adaptive  control  signal  is  different  for  each  actuator,  the  effect 
of  the  MRAS  adaptive  controller  is  a  spatially  varying  gain  on  the  tweeter.  This 
means  that  the  gain  on  each  actuator  varies  as  the  MRAS  controller  adapts  to  the 
deviations  from  the  model.  Because  the  tweeter  commands  are  also  transferred  on  to 
the  woofer  as  seen  in  Figure  2.8,  the  spatially  varying  control  law  is  represented  on 
the  woofer  as  well.  Figure  4.6  shows  the  spatial  variation  of  the  parameter,  so-  The 
parameter,  t0,  also  shows  a  similar  pattern.  Figure  4.7  shows  the  time  history  of  the 
parameter,  s0,  for  one  actuator.  The  figure  suggests  that  the  parameter  has  not  yet 
converged  on  its  final  value  even  after  several  iterations.  This  behavior  is  expected 
since  the  goal  of  the  MRAS  is  to  drive  the  error  to  zero,  not  estimate  parameters  [1]. 
Figure  4.8  shows  the  error  between  the  model  response  and  the  response  of  the  actual 
system  for  a  single  actuator.  After  an  initial  transient,  the  output  of  the  AO  system 
essentially  matches  the  desired  model  output. 

Figure  4.9  compares  the  performance  of  the  MRAS  controller  with  the  baseline 
woofer-tweeter  AO  system.  To  allow  for  a  direct  comparison,  the  same  turbulence 
was  used  in  both  simulations.  To  reduce  the  number  of  figures,  both  the  high  and 


50 


Parameter  Convergence 


Figure  4.7:  Time  history  of  adjustment  parameter  s0- 


MRAS  Modeling  Error 


Figure  4.8:  The  MRAS  error,  y  —  ym,  measures  how  well  the 
true  system  follows  the  desired  model. 
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Comparison  of  MRAS  Filter  and  Baseline 


Figure  4.9:  Comparison  of  fixed-gain  Woofer-Tweeter  AO  and 
MRAS  adaptively  controlled  Woofer- Tweeter. 

low  turbulence  cases  are  included  on  the  same  plot.  In  the  low  turbulence,  the  MRAS 
controller  offers  little  benefit  with  an  increase  in  the  average  Strclil  ratio  of  only 
0.009.  This  increase  in  the  Strclil  ratio  represents  only  a  0.94%  improvement  over  the 
baseline.  The  small  improvement  is  not  surprising  due  to  the  fact  that  the  baseline 
performance  is  already  high  and  leaves  little  room  for  improvement.  However,  it  is 
encouraging  that  the  adaptive  controller  provides  a  benefit  even  in  cases  where  the 
traditional  fixed-gain  controller  offers  good  performance.  This  result  also  validates  the 
adaptive  control  architecture  in  which  only  the  tweeter  is  considered.  On  the  other 
hand,  in  the  presence  of  the  high  turbulence  the  MRAS  controller  provides  significant 
improvement.  The  increase  in  the  average  Strclil  ratio  is  0.073,  or  15.78%.  Also, 
in  both  cases  the  adaptive  controller  provides  the  additional  benefit  of  improving  the 
transient  performance  of  the  AO  system  and  reduces  the  initial  dip  in  the  Strehl  ratio. 

4-3  Lattice  Filter  Simulation  Results 

4-3.1  Actuator- Space.  In  order  to  compare  the  lattice  filter-based  controller 

to  the  baseline,  it  is  necessary  to  determine  both  the  number  of  stages,  or  order,  of 
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Determining  Filter  Order 


Figure  4.10:  The  norm  of  the  forward  error  vector  decreases 
with  filter  order.  The  desired  filter  order  can  be  determined 
by  finding  the  lowest  filter  order  which  still  results  in  a  small 
forward  error  norm. 

the  filter  as  well  as  the  amount  of  training  time  required  prior  to  closing  the  adaptive 
control  loop.  The  order  of  the  filter  was  determined  using  a  method  discussed  by 
Honig  and  Messerschmidt  [13].  In  this  method,  a  vector  Fm  is  created  containing  the 
time  history  of  the  forward  prediction  error  at  the  output  of  each  stage, 

F,n  =  [/(0)m,  /( l)m,  ...,/(n)m]  (4.1) 

The  norm  of  this  forward  error  vector  is  then  computed  for  each  stage  in  the  filter.  As 
the  number  of  stages  increases,  the  norm  of  the  vector  should  decrease.  The  number 
of  stages  can  be  found  by  observing  the  norm  of  the  vector  and  finding  when  the  norm 
begins  to  stop  decreasing.  Figure  4.10  shows  a  plot  of  the  forward  error  vector  norm 
as  a  function  of  filter  order.  Using  these  data  it  was  determined  that  a  filter  order 
of  eight  represented  the  lowest  filter  order  desired.  Using  a  lower  filter  order  would 
result  in  more  error,  while  using  a  higher  filter  order  would  increase  the  computational 
burden  of  the  controller  but  not  significantly  reduce  the  error. 
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Determining  Training  Time 


Figure  4.11:  The  training  time  for  the  lattice  filter  was  deter¬ 
mined  by  plotting  the  mean  Strehl  for  various  training  times. 

In  order  to  determine  the  best  training  time  for  the  lattice  filter,  various  training 
times  were  used.  The  mean  Strehl  ratio  for  each  training  time  was  computed  and 
the  best  training  time  was  then  selected.  When  calculating  the  mean  Strehl  ratio, 
only  the  data  after  closing  the  adaptive  loop  were  considered.  This  ensures  that 
only  the  adaptively  controlled  portion  of  the  data  was  considered  when  measuring 
the  performance.  Figure  4.11  shows  a  plot  of  the  mean  Strehl  ratio  as  a  function 
of  the  number  of  training  samples.  Each  simulation  was  run  for  100  samples.  The 
training  time  refers  to  the  number  of  samples  used  to  train  the  filter  before  closing 
the  adaptive  control  loop.  The  remaining  samples  after  the  training  period,  when  the 
adaptive  controller  was  active,  were  used  to  calculate  the  mean  Strehl  ratio.  From 
this  plot  a  training  time  of  20  samples  was  selected  for  the  lattice  filter  controller. 
When  comparing  the  performance  of  the  lattice  filter  with  that  of  the  fixed-gain  and 
MRAS  controllers,  two  Strehl  ratios  are  computed;  the  first  includes  all  of  the  data 
samples,  and  the  other  includes  only  those  that  occur  after  training.  Excluding  the 
samples  that  occur  prior  to  completing  the  training  focuses  on  the  performance  gained 
only  when  using  the  adaptive  controller.  However,  the  fact  that  the  filter  needs  to 
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Comparison  of  Lattice  and  Baseline 


Figure  4.12:  Comparison  of  lattice  filter  to  baseline. 

be  trained  should  not  be  ignored.  Averaging  over  all  of  the  data  samples,  including 
those  when  the  adaptive  controller  is  not  active,  incorporates  the  training  period  into 
the  performance  metric.  Of  course,  as  the  number  of  data  samples  increases,  the  low 
performance  during  the  training  period  becomes  less  significant. 

With  the  filter  order  specified  and  the  number  of  training  samples  determined, 
the  performance  of  the  lattice  filter  controller  can  be  evaluated.  Figure  4.12  compares 
the  lattice  filter  controller  to  the  fixed-gain  controller.  Prior  to  completion  of  training, 
the  lattice  filter  controller  exhibits  identical  performance  to  the  baseline  fixed-gain 
system  because  without  the  adaptive  control  loop  closed,  only  the  fixed-gain  controller 
is  active.  As  with  the  MRAS  controller,  the  lattice  filter  controller  provides  little 
performance  improvement  when  correcting  for  the  lower  turbulence.  Again,  this  is 
due  to  the  fact  that  there  is  little  room  for  improvement  over  the  fixed-gain  controller. 
The  improvement  in  mean  Strehl  ratio  over  the  entire  data  sample  was  0.025,  or 
2.86%.  When  only  the  active  period  of  the  controller  is  considered  the  mean  Strehl 
ratio  improved  by  0.027  for  a  2.9%  improvement.  However,  in  cases  of  more  severe 
turbulence,  the  performance  improves  dramatically  when  the  adaptive  controller  loop 
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is  enabled.  Immediately  prior  to  closing  the  adaptive  control  loop  the  Strehl  ratio  was 
0.454.  After  closing  the  loop,  the  Strehl  ratio  rapidly  increases  to  0.55.  The  average 
Strehl  ratio  over  all  the  data  improved  by  0.098,  or  21.1%.  After  the  training  period 
the  improvement  was  0.102,  or  21.8%. 

Like  the  MRAS  controller,  the  lattice  filter  controller  also  exhibits  spatially 
varying  gain  parameters.  Both  the  reflection  coefficients  and  the  joint-process  pa¬ 
rameters  vary  spatially.  Figure  4.13  shows  the  spatial  variation  of  the  first  reflection 
coefficient,  K\.  No  direct  comparison  of  the  lattice  filter  controller  parameters  and  the 
MRAS  controller  parameters  can  be  made  due  to  the  vastly  different  structure  of  the 
two  systems.  Still,  in  both  cases  the  spatially  varying  parameters  equate  to  a  spa¬ 
tially  varying  control  law  which  adapts  to  both  the  turbulence  and  any  un-modeled 
effects  of  the  DM  and  WFS.  For  both  the  MRAS  and  the  actuator-space  lattice  filter, 
this  adaptation  occurs  for  each  actuator  based  on  its  time  history.  Figure  4.14  shows 
the  time  history  of  the  parameter,  K\.  for  one  actuator.  After  an  initial  transient, 
the  parameter  converges  toward  a  steady-state  value  and  then  fluctuates  around  that 
value.  This  is  different  than  the  MRAS  controller  where  the  parameters  do  not  reach 
their  steady  state  values.  The  error  signal,  y  —  ym,  measures  how  well  the  true  system 
follows  the  desired  model.  In  the  lattice  filter  controller,  this  signal  is  the  input  to  the 
lattice  filter.  Figure  4.15  shows  a  plot  of  the  the  error  signal  at  one  actuator  location. 

4-3.2  Modal- Space.  Thus  far,  the  lattice  filter  controller  requires  the  es¬ 

timation  of  two  parameters  for  each  stage:  the  reflection  coefficient,  Krn ,  as  well  as 
the  joint-process  weight,  hm.  If  M  is  the  number  of  filter  stages  and  Np  is  the  total 
number  of  tweeter  actuators,  the  total  number  of  parameters,  Np,  estimated  by  the 
filter  is 

Np  =  2  MNt.  (4.2) 

In  this  case,  with  1,245  active  tweeter  actuators  and  8  stages,  the  total  number  of  filter 
parameters  is  19,920.  In  comparison,  the  MRAS  controller  requires  the  estimation 
of  only  two  parameters,  So  and  t0,  for  each  actuator.  This  gives  a  total  of  only 


56 


Spatial  Variaition  of  Parameter  k\ 
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Figure  4.13:  The  spatial  variation  of  parameter  K\  at  each 

tweeter  actuator. 


Convergence  of  Parameter  ki 


Figure  4.14:  The  time  history  of  parameter  K\  for  one  actua¬ 
tor  demonstates  how  the  lattice  filter  parameters  converge  over 
time. 
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Lattice  Modeling  Error 


Figure  4.15:  The  error  signal,  y  —  ym,  measures  how  well  the 
true  system  follows  the  desired  model.  This  signal  is  the  input 
to  the  lattice  filter. 

2,490  parameters,  or  one  eighth  that  of  the  lattice  filter.  The  total  number  of  modes 
computed  in  the  lattice  filter  design  of  Section  3.2.2  is  equal  to  the  number  of  tweeter 
actuators.  Selecting  a  subset  of  the  modes  to  control  reduces  the  number  of  parameters 
and  the  computational  burden.  Figure  4.16  shows  a  plot  of  the  mean  Strehl  ratio  as  a 
function  of  the  number  of  modes.  It  should  be  noted  that  in  order  to  shorten  the  length 
of  the  simulation,  only  50  time  steps  were  used  to  generate  the  plot  in  Figure  4.16. 
Because  the  simulation  time  was  shorter  than  the  other  simulations,  no  comparison 
should  be  made  between  this  plot  and  the  performance  of  the  other  controllers.  The 
purpose  of  the  plot  is  simply  to  identify  the  minimum  number  of  modes  necessary  to 
still  achieve  a  high  level  of  performance.  Based  on  these  data,  800  modes  were  selected 
for  the  modal-space  lattice  filter  controller.  Using  800  modes  reduces  the  number  of 
parameters  by  36%,  from  19,920  to  12,800.  Figure  4.17  compares  the  performance 
of  the  modal-space  lattice  filter  controller  to  the  baseline  case.  As  expected,  the 
improvement  in  the  low  turbulence  scenario  is  minimal  with  an  improvement  of  only 
2%  in  the  mean  Strehl  ratio.  In  the  high  turbulence  scenario  the  improvement  is 
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Performance  vs  Number  of  Modes 


Figure  4.16:  Performance  versus  number  of  modes 


Comparison  of  Modal-Space  Lattice  and  Baseline 


Figure  4.17:  Comparison  of  modal-space  to  baseline. 
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Performance  of  Modal-Space  Lattice  Filter  (155  modes) 


Figure  4.18:  Modal-space  performance  with  155  modes. 

17.1%.  This  performance  is  more  than  the  MRAS  controller,  but  slightly  less  than 
the  actuator-space  controller  using  the  lattice  filter  approach.  The  results  suggest  that 
there  is  a  tradeoff  between  performance  and  computational  burden.  To  illustrate  this 
tradeoff,  the  number  of  modes  in  the  lattice  filter  is  further  reduced  so  the  total  number 
of  lattice  filter  parameters  is  the  same  as  the  number  of  MRAS  parameters.  Setting 
Np  =  2490  and  solving  Equation  4.2  for  Np,  the  required  number  of  modes  is  found  to 
be  155.  Figure  4.18  compares  the  performance  of  the  modal-space  controller  with  the 
reduced  number  of  modes.  While  the  performance  still  exceeds  that  of  the  baseline 
system,  the  improvement  is  substantially  less.  The  mean  Strehl  ratio  improved  by  only 
0.67%  and  4.98%  for  the  weak  and  strong  turbulent  cases,  respectively.  This  represents 
roughly  one  third  the  improvement  seen  by  the  other  controllers,  showing  that  the 
lattice  filter  performance  degrades  when  fewer  modes  are  used.  It  also  shows  that 
when  the  lattice  filter  uses  the  same  number  of  parameters  as  the  MRAS  controller, 
much  of  the  performance  improvement  is  lost  and  the  lattice  filter  controller  no  longer 
outperforms  the  MRAS  controller.  The  degradation  in  the  lattice  filter  controller 
performance  is  due  to  the  fact  that  with  fewer  modes  more  information  is  being 
discarded  by  the  adaptive  controller.  This  implies  that  if  the  controller  has  a  limited 
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computational  ability,  the  MRAS  controller  can  result  in  better  performance  than 
the  lattice  filter  controller,  clue  to  the  fact  that  the  MRAS  controller  does  not  discard 
information. 

With  the  simulation  of  the  different  controllers,  it  is  clear  that  in  all  of  the 
cases  examined  that  adaptive  control  improves  the  performance  of  woofer-tweeter 
AO  systems.  Chapter  V  provides  a  brief  summary  of  the  results  presented  here  and 
gives  recommendations  for  future  research  on  the  topic  of  adaptive  control  for  woofer- 
tweeter  AO. 
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V.  Conclusions  and  Recommendations 


This  chapter  summarizes  the  research  effort  in  adaptive  control  of  woofer-tweeter 
AO  and  the  simulation  results  presented  in  Chapter  IV.  The  results  are  used  to 
provide  recommendations  for  further  research.  These  recommendations  aim  to  further 
develop  the  control  algorithms  presented  and  extend  them  to  other  research  areas  of 
interest  to  the  AO  community. 

5.1  Summary  and  Conclusion 

The  goal  of  this  research  was  to  develop  an  adaptive  control  algorithm  for  a 
woofer-tweeter  AO  system  and  to  verify  its  performance  with  wave-optics  simulations. 
The  purpose  of  the  adaptive  controller  was  to  improve  AO  performance,  particularly 
in  the  presence  of  high- variance,  fast-moving  turbulence. 

The  first  step  in  the  design  process  was  to  determine  a  general  control  archi¬ 
tecture  for  implementing  the  adaptive  controller.  Based  on  recent  work  by  Gibson  et 
al.  [8,9,14-17,19],  a  control  architecture  where  the  adaptive  controller  provides  addi¬ 
tional  input  to  the  standard  fixed-gain  controller  was  used.  This  control  architecture 
was  adapted  to  a  woofer-tweeter  AO  system  by  applying  the  adaptive  controls  only  to 
the  tweeter  and  then  transferring  the  average  tweeter  commands  to  the  woofer.  Using 
this  control  architecture,  three  different  adaptive  controllers  were  developed:  a  model 
reference  adaptive  system  controller,  a  lattice  filter  controller  operating  in  actnator- 
space,  and  a  lattice  filter  controller  operating  in  modal-space.  The  comparatively 
simple  MRAS  controller  was  used  to  validate  the  general  control  architecture.  The 
lattice  filter  controllers  were  developed  to  overcome  the  potential  instabilities  inherent 
in  MRAS  controllers.  A  modal  decomposition  of  the  DM  commands  was  introduced 
to  both  reduce  the  number  of  filter  parameters  and  allow  the  filter  to  adapt  to  both 
spatial  and  temporal  variations  in  the  system. 

After  designing  the  controllers,  they  were  simulated  in  two  different  scenarios, 
one  with  high  and  one  with  low  turbulence,  to  measure  their  performance.  A  tracli- 


62 


Table  5.1:  Summary  of  Results:  For  each  turbulence  scenario  and  each 

controller  type,  the  number  of  adaptation  parameters,  average  Strehl  ratio, 
and  percentage  change  from  the  baseline  Strehl  ratio  are  shown. 


Scenario 

Controller 

NP 

Strehl 

%  Change 

Low  Turb 

Fixed-Gain 

0 

0.9063 

- 

MRAS 

2,490 

0.9148 

0.94% 

Lattice  (Actuator-Space) 

19,920 

0.9370 

2.83% 

Lattice  (800  Modes) 

12,800 

0.9245 

2.00% 

Lattice  (155  Modes) 

2,480 

0.9124 

0.67% 

High  Turb 

Fixed-Gain 

0 

0.4625 

- 

MRAS 

2,490 

0.5355 

15.78% 

Lattice  (Actuator-Space) 

19,920 

0.5602 

21.13% 

Lattice  (800  Modes) 

12,800 

0.5416 

17.10% 

Lattice  (155  Modes) 

2,480 

0.4855 

4.98% 

tional  fixed-gain  controller  was  used  as  the  baseline  case  for  each  scenario.  Table  5.1 
summarizes  the  simulation  results  for  each  of  the  adaptive  controllers  used. 

In  the  low  turbulence  scenario,  none  of  the  adaptive  controllers  resulted  in  a 
substantial  performance  increase.  The  actuator-space  lattice  filter  controller,  which 
resulted  in  the  most  improvement,  only  increased  the  average  Strehl  ratio  by  2.8%. 
These  small  gains  in  low  turbulence  are  clue  to  the  fact  that  the  overall  AO  perfor¬ 
mance  is  already  high  and  does  not  leave  much  room  for  improvement.  However,  in 
all  cases  the  adaptive  controllers  improved  the  performance  of  the  woofer-tweeter  AO 
system. 

The  simulation  results  in  the  high  turbulence  scenario  demonstrate  the  perfor¬ 
mance  differences  between  the  various  adaptive  controllers.  Again,  the  actuator-space 
lattice  filter  controller  showed  the  most  improvement,  with  a  21%  increase  in  the  av¬ 
erage  Strehl  ratio.  This  performance  comes  with  a  computational  burden  of  nearly 
20,000  filter  parameters.  If  the  number  of  parameters  is  reduced  by  roughly  30%, 
through  the  use  of  a  modal  decomposition,  the  average  Strehl  ratio  still  increases  by 
17%.  If,  however,  the  number  of  lattice  filter  parameters  is  further  reduced  to  be 
equal  to  the  number  of  MRAS  parameters,  the  modal-space  lattice  filter  controller 
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only  increases  the  average  Strehl  ratio  by  5%,  while  the  MRAS  controller  acheives  a 
17%  improvement  with  the  same  number  of  parameters. 

In  summary,  the  actuator-space  lattice  filter  adaptive  controller  results  in  the 
most  improvement  in  AO  performance,  in  both  low  and  high  turbulence.  This  per¬ 
formance  comes  at  a  relatively  large  computational  cost  as  measured  by  the  number 
of  filter  parameters  required.  If  a  similar  number  of  parameters  are  used  in  both  the 
lattice  filter  and  the  MRAS  controllers,  the  MRAS  controller  results  in  better  perfor¬ 
mance.  However,  the  MRAS  controller  does  not  guarantee  stability.  A  poorly  chosen 
step  size  parameter,  7,  could  possibly  lead  to  instability.  The  lattice  filter,  on  the 
other  hand,  with  a  FIR  implementation,  guarantees  BIBO  stability  and  has  no  arbi¬ 
trary  parameters  that  can  affect  stability.  For  these  reasons,  the  modal-space  lattice 
filter  controller  with  800  modes  is  the  best  controller  for  the  scenarios  presented.  It 
offers  good  performance,  guaranteed  stability,  and  a  lower  number  of  filter  parame¬ 
ters.  In  practice,  each  AO  system  is  unique  and  none  of  the  presented  controllers  can 
be  considered  superior  to  the  others  in  all  cases.  In  fact,  if  only  a  limited  amount 
of  computational  resources  are  available,  the  MRAS  controller  would  be  a  reasonable 
choice. 

5. 2  Recommendations 

Several  alternatives  exist  for  implementing  or  improving  upon  the  adaptive  con¬ 
trol  algorithms  presented.  This  section  discusses  some  of  these  areas  of  possible  future 
research. 

The  results  presented  in  this  research  were  all  based  solely  on  computer  simu¬ 
lations.  Perhaps,  the  most  logical  extension  of  this  research  is  to  implement  and  test 
the  controller  designs  using  real  hardware.  In  any  simulation  it  is  impossible  to  fully 
capture  the  intricacies  of  the  real  world.  With  real  devices,  though,  all  un-modeled 
effects  are  unavoidably  included  and  it  is  possible  to  obtain  a  true  measure  of  the 
system  performance.  Prior  to  implementing  the  adaptive  controller  with  real  hard¬ 
ware  it  is  possible  to  predict  the  average  improvement  and  standard  deviation  of  a 
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given  controller  using  Monte  Carlo  simulations  using  different  atmospheric  turbulence 
seeds. 

Another  possible  area  of  future  research  would  be  to  implement  the  controller 
with  more  complicated  systems.  As  mentioned  briefly  in  Chapter  II,  some  recent 
work  has  been  done  using  spatially  varying  fixed-gain  controllers  [21].  If  a  spatially- 
varying  gain  were  used  on  the  fixed-gain  controller,  these  same  gains  would  need  to 
be  incorporated  into  the  system  model  used  in  the  adaptive  controllers.  While  this 
change  may  be  fairly  straightforward  in  the  controllers  which  operate  in  actuator 
space,  research  would  have  to  be  done  on  how  to  map  the  spatially  varying  gain  into 
the  modal-space  model.  It  is  conceivable  that  such  research  could  yield  an  adaptive 
control  law  where  the  gains  of  the  model  system  vary  by  mode  instead  of  modelling 
each  mode  identically. 

Another  possibility  would  be  to  investigate  the  contribution  of  each  mode  over 
time.  If  certain  modes  are  found  to  be  more  significant  sources  of  error,  it  might  be 
possible  to  design  an  adaptive  controller  which  applies  itself  only  to  the  significant 
modes.  Such  a  controller  might  be  able  to  maintain  a  high  level  of  performance  with 
fewer  adaptive  parameters  by  dynamically  selecting  which  modes  to  control. 

The  last  recommendation  primarily  relates  to  fixed-gain  woofer- tweeter  AO.  It 
is  possible  that  using  an  offloading  technique  other  than  simply  downsampling  the 
tweeter  command,  might  change  the  performance  of  the  woofer-tweeter  AO  system. 
An  example  of  how  this  might  be  done  would  be  to  generate  the  woofer  commands 
based  on  an  weighted  average  of  the  adjacent  tweeter  commands  by  creating  a  filtered 
version  of  Ht2W-  Since  the  offload  matrix  did  not  factor  into  the  design  of  the  adaptive 
controllers,  this  change  would  not  effect  the  design  of  the  adaptive  controller. 

These  recommendations  represent  only  a  few  of  the  areas  where  future  research 
can  be  conducted.  There  are  likely  several  possible  areas  of  research  that  have  not 
been  considered  here.  Given  the  demonstrated  benefits  of  adaptive  control  to  woofer- 
tweeter  AO  systems,  further  research  in  the  area  is  encouraged. 
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